CHAPTER 



Oscillations 



Almost any system that is displaced from a position of stable equilibrium exhibits 
oscillations . If the displacement is small , the oscillations are almost always of the type 
called simple harmonic. Oscillations, and particularly simple harmonic oscillations, 
are therefore extremely widespread. They are also extremely useful. For example, 
all good clocks depend on an oscillator to regulate their time keeping: The first 
reliable clocks used a pendulum; the first accurate watches (historically crucial in 
navigation) used an oscillating balance wheel; modern watches use the oscillations 
of a quartz crystal; and today’s most accurate clocks, such as the atomic clock at 
the National Institute for Standards and Technology in Boulder, Colorado, use the 
oscillations of an atom. In this chapter, we shall explore the physics and mathematics 
of oscillations. I shall begin with simple harmonic oscillations and then go on to 
damped oscillations (oscillations that die out because of resistive forces) and driven 
oscillations (oscillations that are maintained by an outside driving force, as in all 
clocks). The last three sections of this chapter describe the use of Fourier series in 
finding the motion of an oscillator driven by an arbitrary periodic driving force. 


5.1 Hooke’s Law 


As you are certainly aware, a mass on the end of a spring that obeys Hooke’s law 
executes oscillations of the type that we call simple harmonic. Before we review the 
proof of this claim, let us first ask why Hooke’s law is so important and appears so 
frequently. Hooke’s law asserts that the force exerted by a spring has the form (for 
now we’ll restrict ourselves to a spring confined to the x axis) 

F x (x) = —kx (5.1) 

where x is the displacement of the spring from its equilibrium length and k is a positive 
number called the force constant. That k is positive means that the equilibrium at x = 0 
is stable: When x = 0 there is no force, when x > 0 (displacement to the right) the 161 
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force is negative (back to the left), and when x < 0 (displacement to the left) the 
force is positive (back to the right); either way, the force is a restoring force, and the 
equilibrium is stable. (If k were negative, the force would be away from the origin, and 
the equilibrium would be unstable, in which case we do not expect to see oscillations.) 
An exactly equivalent way to state Hooke’s law is that the potential energy is 

U(x) = \kx 2 . 


Consider now an arbitrary conservative one-dimensional system which is specified 
by a coordinate x and has potential energy U(x). Suppose that the system has a stable 
equilibrium position x = x 0 , which we may as well take to be the origin (x 0 = 0). 
Now consider the behavior of U (x) in the vicinity of the equilibrium position. Since 
any reasonable function can be expanded in a Taylor series, we can safely write 

U(x) = U(0) + U'(0)x + \U"{ 0)x 2 + ■■■. (5.2) 

As long as x remains small, the first three terms in this series should be a good 
approximation. The first term is a constant, and, since we can always subtract a 
constant from U (x) without affecting any physics, we may as well redefine U (0) 
to be zero. Because x = 0 is an equilibrium point, U'( 0) = 0 and the second term 
in the series (5.2) is automatically zero. Because the equilibrium is stable, U'\ 0) is 
positive. Renaming U"( 0) as k , we conclude that for small displacements it is always 
a good approximation to take 1 

U{x) = \kx 2 . (5.3) 

That is, for sufficiently small displacements from stable equilibrium, Hooke’s law is 
always valid. Notice that if U"( 0) were negative, then k would also be negative, and the 
equilibrium would be unstable — a case we’re not interested in just now. Hooke’s law 
in the form (5.3) crops up in many situations, although it is certainly not necessary that 
the coordinate be the rectangular coordinate x, as the following example illustrates. 
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example 5.1 The Cube Balanced on a Cylinder 

Consider again the cube of Example 4.7 (page 130) and show that for small 
j angles 6 the potential energy takes the Hooke’s-law form U(9) = jkO 2 . 

We saw in that example that 

I 

l 

j U(6) = mg[(r + b) cos 6 + r6 sin#]. 

I 

| 

I 

| If 0 is small we can make the approximations cos 0 ~ 1 — 0 2 /2 and sin 0 ~ 0, 
so that 

] 

I U(0) ~ mg[(r + b)( 1 — 1 0 2 ) + rO 2 ] = mg(r + b) + \mg(r — b)6 2 , 


1 The only exception is if U"( 0) happens to be zero, but I shall not worry about this exceptional 
case here. 
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which, apart from the uninteresting constant, has the form \k0 2 with “spring 
constant” k = mg(r — b). Notice that the equilibrium is stable (k positive) only 
when r > b, a condition we had already found in Example 4.7. 



Figure 5.1 A mass m with potential energy U(x) = \kx 2 and total 
energy E oscillates between the two turning points at x = ±A, where 
U(x) = E and the kinetic energy is zero. 


As discussed in Section 4.6, the general features of the motion of any one¬ 
dimensional system can be understood from a graph of U(x ) against x. For the 
Hooke’s-law potential energy (5.3), this graph is a parabola, as shown in Figure 
5.1. If a mass m has potential energy of this form and has any total energy E > 0, 
it is trapped and oscillates between the two turning points where U(x) = E, so that 
the kinetic energy is zero and the mass is instantaneously at rest. Because U(x ) is 
symmetric about x = 0, the two turning points are equidistant on opposite sides of 
the origin and are traditionally denoted x — ±A, where A is called the amplitude of 
the oscillations. 


5.2 Simple Harmonic Motion 


We are now ready to examine the equation of motion (that is, Newton’s second 
law) for a mass m that is displaced from a position of stable equilibrium. To be 
definite, let us consider a cart on a frictionless track attached to a fixed spring as 
sketched in Figure 5.2. We have seen that we can approximate the potential energy 
by (5.3) or, equivalently, the force by F x {x) = —kx. Thus the equation of motion is 
mx = F x — —kx or 

x = — —x = —co 2 x (5.4) 

m 

where I have introduced the constant 



which we shall see is the angular frequency with which the cart will oscillate. Although 
we have arrived at Equation (5.4) in the context of a cart on a spring moving along 
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Figure 5.2 A cart of mass m oscillating on the end of a spring. 


the x axis, we shall see eventually that it applies to many different oscillating systems 
in many different coordinate systems. For example, we have already seen in Equation 
(1.55) that the angle 0 that gives the position of a pendulum (or a skateboard in a 
half-pipe) is governed by the same equation, 0 = —&> 2 0, at least for small values 
of 0. In this section I am going to review the properties of the solutions of (5.4). 
Unfortunately, there are many different ways to write the same solution, all of which 
have their advantages, and you should be comfortable with them all. 


The Exponential Solutions 

Equation (5.4) is a second-order, linear, homogeneous differential equation 2 and so 
has two independent solutions. These two independent solutions can be chosen in 
several different ways, but perhaps the most convenient is this: 

x(t) = e l(ot and x(t) = e ~ lC0t . 

As you can easily check, both of these functions do satisfy (5.4). Further, any constant 
multiple of either solution is also a solution, and likewise any sum of such multiples. 
Thus the function 


x(t) = C x e i(0t + C 2 e~ icot (5.5) 

is also a solution for any two constants C l and C 2 . (That any linear combination 
of solutions like this is itself a solution is called the superposition principle and 
plays a crucial role in many branches of physics.) Since this solution (5.5) contains 
two arbitrary constants, it is the general solution of our second-order equation (5.4). 3 
Therefore, any solution can be expressed in the form (5.5) by suitable choice of the 
coefficients C x and C 2 . 


2 Linear because it contains no higher powers of x or its derivatives than the first power, and 
homogeneous because every term is a first power (that is, there is no term independent of x and its 
derivatives). 

3 Recall the result, discussed below Equation (1.56), that the general solution of a second-order 
differential equation contains precisely two arbitrary constants. 
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The Sine and Cosine Solutions 

The exponential functions in (5.5) are so convenient to handle that (5.5) is often the 
best form of the solution. Nevertheless, this form does have one disadvantage. We 
know, of course, that x(t) is real, whereas the two exponentials in (5.5) are complex. 
This means the coefficients C x and C 2 must be chosen carefully to ensure that x(t) 
itself is real. I shall return to this point shortly, but first I shall rewrite (5.5) in another 
useful way. From Euler’s formula (2.76) we know that the two exponentials in (5.5) 
can be written as 


e ±lC0t = cos (cot) d= i sin (cot). 

Substituting into (5.5) and regrouping we find that 

x(t) = (Ci + C 2 ) cos (cot) + i(Ci — C 2 ) sin(cot) 

= B x cos (cot) + B 2 sin(&>0 (5.6) 

where B x and B 2 are simply new names for the coefficients in the previous line, 

B x = Ci + C 2 and B 2 = i(C x - C 2 ). (5.7) 

The form (5.6) can be taken as the definition of simple harmonic motion (or SHM): 
Any motion that is a combination of a sine and cosine of this form is called simple 
harmonic. Because the functions cos (cot) and sin (cut) are real, the requirement that 
x(t ) be real means simply that the coefficients B ] and B 2 must be real. 

We can easily identify the coefficients B x and B 2 in terms of the initial conditions 
of the problem. Clearly at t = 0, (5.6) implies that jt(0) = B x . That is, B x is just the 
initial position v(0) = x 0 . Similarly, by differentiating (5.6), we can identify coB 2 as 
the initial velocity v 0 . 

If I start the oscillations by pulling the cart aside to x = x 0 and releasing it from 
rest (v 0 = 0), then B 2 — 0 in (5.6) and only the cosine term survives, so that 

x(t) = x 0 cos(cot). (5.8) 

If I launch the cart from the origin (x 0 — 0) by giving it a kick at t = 0, only the sine 
term survives, and 


v 

x(t) = — sin(ct>0- 
co 

These two simple cases are illustrated in Figure 5.3. Notice that both solutions, like 
the general solution (5.6), are periodic because both the sine and cosine are. Since the 
argument of both sine and cosine is cot , the function x(t ) repeats itself after the time 
t for which cox — 2n. That is, the period is 



r = 


co 


(5.9) 
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Figure 5.3 (a) Oscillations in which the cart is released from x Q at t = 0 follow 

a cosine curve, (b) If the cart is kicked from the origin at t = 0, the oscillations 
follow a sine curve with initial slope v Q . In either case the period of the oscilla¬ 
tions is r = 2jt/cjo = In^/m/ k and is the same whatever the values of x 0 or v 0 . 


The Phase-Shifted Cosine Solution 

The general solution (5.6) is harder to visualize than the two special cases of Fig¬ 
ure 5.3, and it can usefully be rewritten as follows: First, we define yet another 
constant 


A = 



(5.10) 


Notice that A is the hypotenuse to a right triangle whose other two sides are B x and 
B 2 .1 have indicated this in Figure 5.4, where I have also defined S as the lower angle 
of that triangle. We can now rewrite (5.6) as 


x(t ) = A 


B 1 B 2 

— cos(cot) H-sin(<ztf) 

A A 


= A[ cos 8 cos (cot) + sin 8 sin(&tf)] 
= A cos (cot — 8). 


(5.11) 


From this form it is clear that the cart is oscillating with amplitude A, but instead of 
being a simple cosine as in (5.8), it is a cosine which is shifted in phase: When t = 0 
the argument of the cosine is —8, and the oscillations lag behind the simple cosine by 
the phase shift 8. We have derived the result (5.11) from Newton’s second law, but, as 
so often happens, one can derive the same result in more than one way. In particular, 



Figure 5.4 The constants A and 8 are defined in terms of B x and 
B 2 as shown. 
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(5.11) can also be derived using the energy approach discussed in Section 4.6. (See 
Problem 4.28.) 

Solution as the Real Part of a Complex Exponential 

There is still another useful way to write our solution, in terms of the complex 
exponentials of (5.5). The coefficients C 1 and C 2 there are related to the coefficients 
B l and B 2 of the sine-cosine form by Equation (5.7), which we can solve to give 

C l = \{B l -iB 2 ) and C 2 = \{B X +iB 2 ). (5.12) 

Since B x and B 2 are real, this shows that both C x and C 2 are generally complex and 
that C 2 is the complex conjugate of C h 

c 2 = c*. 

(Recall that for any complex number z = x + iy , the complex conjugate z* is defined 
as 4 z* = x — iy.) Thus the solution (5.5) can be written as 

x(t) = C x e i(0t + C*e~ ieot (5.13) 

where the whole second term on the right is just the complex conjugate of the first. 
(See Problem 5.35 if this isn’t clear to you.) Now, for any complex number z = x + iy , 

z + z* = (x 4- iy) + (jc — iy) = 2x = 2 Re z 

where Rez denotes the real part of z (namely x). Thus (5.13) can be written as 

x(t) = 2Re<V'^. 

If we define a final constant C = 2C b we see from Equation (5.12) and Figure 5.4 
that 


C = B l -iB 2 = Ae~ iS 


(5.14) 


and 


x(t) = Re Ce Uot = Re Ae i(a)t ~ S) . 

This beautiful result is illustrated in Figure 5.5. The complex number Ac' 1 "’ 1 moves 
counterclockwise with angular velocity co around a circle of radius A. Its real part 
[namely x(t)] is the projection of the complex number onto the real axis. While the 
complex number goes around the circle, this projection oscillates back and forth on the 
x axis, with angular frequency co and amplitude A. Specifically, x(t) = A cos (cot — 8), 
in agreement with (5.11). 


4 While most physicists use the notation z*, mathematicians almost always use z for the complex 
conjugate of z. 
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Figure 5.5 The position x(0 of the oscillating cart is the real 
part of the complex number Ae l . As the latter moves around 

the circle of radius A, the former oscillates back and forth on the 
x axis with amplitude A. 


example 5.2 A Bottle in a Bucket 

A bottle is floating upright in a large bucket of water as shown in Figure 5.6. In 
equilibrium it is submerged to a depth d Q below the surface of the water. Show 
that if it is pushed down to a depth d and released, it will execute harmonic 
motion, and find the frequency of its oscillations. If d Q = 20 cm, what is the 
period of the oscillations? 

The two forces on the bottle are its weight mg downward and the upward 
j buoyant force QgAd , where q is the density of water and A is the cross-sectional 
I area of the bottle. (Remember that Archidemes’ principle says that the buoyant 
force is Qg times the volume submerged, which is just Ad.) The equilibrium 
depth d 0 is determined by the condition 

j mg = QgAd 0 . (5.15) 



Figure 5.6 The bottle shown has been loaded with sand so that 
it floats upright in a bucket of water. Its equilibrium depth is 
d — d 0 . 
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Suppose now the bottle is at a depth d = d 0 + x. (This defines x as the distance 
from equilibrium , always the best coordinate to use.) Newton’s second law now 
reads 

mi = mg — QgA(d 0 + x). 


By (5.15), the first and second terms on the right cancel, and we’re left with 
x = —QgAx/m. But again by (5.15), QgA/m = g/d 0 , so the equation of motion 
becomes 


x = — — x, 


which is exactly the equation for simple harmonic motion. We conclude that 
the bottle moves up and down in SHM with angular frequency co — gg/d o- A 
remarkable feature of this result is that the frequency of oscillations is indepen¬ 
dent of m, q, and A; also, the frequency is the same as that of a simple pendulum 
of length / = d 0 . If d 0 — 20 cm, then the period is 


lit 


d r 


t — — =2n / — = 2tt 


co 


g 


0.20 m 
9.8 m/s 2 


= 0.9 sec. 


Try this experiment yourself! But be aware that the details of the flow of water 
around the bottle complicate the situation considerably. The calculation here is 
a very simplified version of the truth. 


Energy Considerations 

To conclude this section on simple harmonic motion, let us consider briefly the energy 
of the oscillator (the cart on a spring or whatever it is) as it oscillates back and forth. 
Since x(t) = A cos (cot — 8 ), the potential energy is just 

U = \kx 2 = \kA 2 cos 2 (cot — 8). 

Differentiating x(t ) to give the velocity, we find for the kinetic energy 

T = jmi 2 = \mo? A 1 sin 2 (&tf — 8) 

= \kA 2 sin 2 (o;? — <5) 

where the second line results from replacing co 2 by k/m. We see that both U and T 
oscillate between 0 and \kA 2 , with their oscillations perfectly out of step — when U 
is maximum T is zero and vice versa. In particular, since sin 2 6 + cos 2 0 = 1, the total 
energy is constant, 


E = T + U = \kA 2 , 


(5.16) 


as it has to be for any conservative force. 
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5.3 Two-Dimensional Oscillators 


In two or three dimensions, the possibilities for oscillation are considerably richer 
than in one dimension. The simplest possibility is the so-called isotropic harmonic 
oscillator, for which the restoring force is proportional to the displacement from 
equilibrium, with the same constant of proportionality in all directions: 

F = —kr. (5.17) 


That is, F x = —kx, F y — —ky (and F z = — kz in three dimensions), all with the same 
constant k. This force is a central force directed toward the equilibrium position, which 
we may as well take to be the origin, as sketched in Figure 5.7(a). Figure 5.7(b) shows 
an arrangement of four identical springs that would produce a force of this form; it is 
easy to see that if the mass at the center is moved away from its equilibrium position it 
will experience a net inward force, and it is not too hard to show (Problem 5.19) that 
this inward force has the form (5.17) for small displacements r . 5 Another example of 
a two-dimensional isotropic oscillator is (at least approximately) a ball bearing rolling 
near the bottom of a large spherical bowl. Two important three-dimensional examples 
are an atom vibrating near its equilibrium in a symmetric crystal, and a proton (or 
neutron) as it moves inside a nucleus. 

Let us consider a particle that is subject to this type of force and suppose, for 
simplicity, that it is confined to two dimensions. The equation of motion, r = F/ra, 
separates into two independent equations: 


2 1 

X = —cox 
2 

y = —coy 


(5.18) 


where I have introduced the familiar angular frequency oo = k/m (which is the same 
in both x and y equations because the same is true of the force constants). Each of 



Figure 5.7 (a) A restoring force that is proportional to r defines 

the isotropic harmonic oscillator, (b) The mass at the center of this 
arrangement of springs would experience a net force of the form 
F = — kr as it moves in the plane of the four springs. 


5 It is perhaps worth pointing out that one does not get a force of the form (5.17) by simply 
attaching a mass to a spring whose other end is anchored to the origin. 
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these two equations has exactly the form of the one-dimensional equation discussed 
in the last section, and the solutions are [as in (5.11)] 

x(t) = A x cos(cot — 8 X ) 
y(t ) = A y cos(<jot — 8 y ) 

where the four constants A x , A y , 8 X , and 8 y are determined by the initial conditions 
of the problem. By redefining the origin of time, we can dispose of the phase shift 8 X , 
but, in general, we cannot also dispose of the corresponding phase in the y solution. 
Thus the simplest form for the general solution is 

x(t) = A x cos(cot ) 
y(t) = A y cos(a)t — 8) 

where 8 = 8 y — 8 X is the relative phase of the y and x oscillations. (See Problem 5.15.) 

The behavior of the solution (5.20) depends on the values of the three constants 
A x , A y , and 8. If either A x or A y is zero, the particle executes simple harmonic motion 
along one of the axes. (The ball bearing in the bowl rolls back and forth through the 
origin, moving in the x direction only or the y direction only.) If neither A x nor A y 
is zero, the motion depends critically on the relative phase 8. If 8 = 0, then x(t) and 
y(t) rise and fall in step, and the point (jc, y) moves back and forth on the slanting 
line that joins (A x , A y ) to (—A x , —A y ), as shown in Figure 5.8(a). If 8 = jt/ 2, then 
jc and y oscillate out of step, with x at an extreme when y is zero, and vice versa; the 
point (jc, y) describes an ellipse with semimajor and semiminor axes A x and A y , as in 
Figure 5.8(b). For other values of 8, the point (jc, y) moves around a slanting ellipse, 
as shown for the case 8 = tt/4 in Figure 5.8(c). (For a proof that the path really is an 
ellipse, see Problem 8.11.) 

In the anisotropic oscillator, the components of the restoring force are propor¬ 
tional to the components of the displacement, but with different constants of propor¬ 
tionality: 

F x = -k x x, F y = -k y y, and F z = -k z z. (5.21) 


(5.20) 


(5.19) 



(a) 8 = 0 (b) S = tt/2 (c) S = n/4 


Figure 5.8 Motion of a two-dimensional isotropic oscillator as given by (5.20). (a) If 
5 = 0, then x and y execute simple harmonic motion in step, and the point (x, y) moves 
back and forth along a slanting line as shown, (b) If 8 = 7t/2, then (x, y) moves around 
an ellipse with axes along the x and y axes, (c) In general (for example, 8 = tt/4), the 
point (x, y) moves around a slanted ellipse as shown. 
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An example of such a force is the force felt by an atom displaced from its equilibrium 
position in a crystal of low symmetry, where it experiences different force constants 
along the different axes. For simplicity, I shall again consider a particle in two 
dimensions, for which Newton’s second law separates into two separate equations 
just as in (5.18): 


** 2 1 

x = —cox 

JC 

y = . 


(5.22) 


The only difference between this and (5.18) is that there are now different frequencies 
for the different axes, co x = y/k x /m and so on. The solution of these two equations is 
just like (5.20): 


x(t) = A x cos(co x t) ^ 5 23 ^ 

y(t) = A y cos(co y t - 8). 

Because of the two different frequencies, there is a much richer variety of possible 
motions. If co x /co y is a rational number, it is fairly easy to see (Problem 5.17) that the 
motion is periodic, and the resulting path is called a Lissajous figure (after the French 
physicist Jules Lissajous, 1822-1880). For example, Figure 5.9(a) shows an orbit of 
a particle for which co x /(jo y — 2 and the x motion repeats itself twice as often as the 
y motion. In the case shown, the result is a figure eight. If co x /co y is irrational, the 
motion is more complicated and never repeats itself. This case is illustrated, with 
(o x /co y = V2, in Figure 5.9(b). This kind of motion is called quasiperiodic: The 
motion of the separate coordinates x and y is periodic, but because the two periods 
are incompatible, the motion of r = (x, y) is not. 



Figure 5.9 (a) One possible path for an anisotropic oscillator with co x = 2 

and co y = 1. You can see that x goes back and forth twice in the time that y 
does so once, and the motion then repeats itself exactly, (b) A path for the 
case co x = y/2 and co y = 1 from t = 0 to t = 24. In this case the path never 
repeats itself, although, if we wait long enough, it will come arbitrarily close 
to any point in the rectangle bounded by x = ±A X and y = ±A y . 
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5.4 Damped Oscillations 


I shall now return to the one-dimensional oscillator, and take up the possibility that 
there are resistive forces that will damp the oscillations. There are several possibilities 
for the resistive force. Ordinary sliding friction is approximately constant in magni¬ 
tude, but always directed opposite to the velocity. The resistance offered by a fluid, 
such as air or water, depends on the velocity in a complicated way. However, as we 
saw in Chapter 2, it is sometimes a reasonable approximation to assume that the re¬ 
sistive force is proportional to v or (under different circumstances) to v 2 . Here I shall 
assume that the resistive force is proportional to v ; specifically, f = — b\. One of my 
main reasons is that this case leads to an especially simple equation to solve, and the 
equation is itself a very important equation that appears in several other contexts and 
is therefore well worth studying. 6 

Consider, then, an object in one dimension, such as a cart attached to a spring, that 
is subject to a Hooke’s law force, — kx, and a resistive force, —bx. The net force on 
the object is — bx — kx , and Newton’s second law reads (if I move the two force terms 
over to the left side) 


mx + bx + kx = 0. (5.24) 

One of the beautiful things about physics is the way the same mathematical 
equation can arise in totally different physical contexts, so that our understanding of 
the equation in one situation carries over immediately to the other. Before we set about 
solving Equation (5.24), I would like to show how the same equation appears in the 
study of LRC circuits. An LRC circuit is a circuit containing an inductor (inductance 
L), a capacitor (capacitance C), and a resistor (resistance R ), as sketched in Figure 
5.10.1 have chosen the positive direction for the current to be counterclockwise, and 
the charge q(t) to be the charge on the left-hand plate of the capacitor [with —q(t) on 
the right], so that I (t) = q (t ). If we follow around the circuit in the positive direction, 
the electric potential drops by LI = Lq across the inductor, by RI = Rq across the 
resistor, and by q/C across the capacitor. Applying Kirchoff’s second rule for circuits, 
we conclude that 


Lq + Rq + — q = 0. (5.25) 

C 

This has exactly the form of Equation (5.24) for the damped oscillator, and anything 
that we learn about the equation for the oscillator will be immediately applicable to 
the LRC circuit. Notice that the inductance L of the electric circuit plays the role of 
the mass of the oscillator, the resistance term Rq corresponds to the resistive force, 
and 1/C to the spring constant k. 


6 You should be aware, however, that although the case I am considering — that the resistive 
force f is linear in v — is very important, it is nevertheless a very special case. I shall describe some 
of the startling complications that can occur when f is not linear in v in Chapter 12 on nonlinear 
mechanics and chaos. 
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Figure 5.10 An LRC circuit. 


Let us now return to mechanics and the differential equation (5.24). To solve this 
equation it is convenient to divide by m and then introduce two other constants. I shall 
rename the constant b/m as 2)6, 

- = 2p. (5.26) 

m 

This parameter /?, which can be called the damping constant, is simply a convenient 
way to characterize the strength of the damping force — as with b , large f corresponds 
to a large damping force and conversely. I shall rename the constant k/m as co 7 , 
that is, 

(5.27) 

Notice that co 0 is precisely what I was calling co in the previous two sections. I have 
added the subscript because, once we admit resistive forces, various other frequencies 
become important. From now on, I shall use the notation co Q to denote the system’s 
natural frequency, the frequency at which it would oscillate if there were no resistive 
force present , as given by (5.27). With these notations, the equation of motion (5.24) 
for the damped oscillator becomes 



X+20X + 


(5,2:1 


Notice that both of the parameters f> and co Q have the dimensions of inverse time, that 
is, frequency. 

Equation (5.28) is another second-order, linear, homogeneous equation [the last 
was (5.4)]. Therefore, if by any means we can spot two independent 7 solutions, x x (t) 
and x 2 (t) say, then any solution must have the form Cpcfif) + C 2 x 2 (t). What this 


7 It is about time I gave you a definition of “independent.” In general this is a little complicated, 
but for two functions it is easy: Two functions are independent if neither is a constant multiple of 
the other. Thus the two functions sin(x) and cos(x) are independent; likewise the two functions x 
and x 2 \ but the two functions .x and 3x are not. 
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means is that we are free to play a game of inspired guessing to find ourselves two 
independent solutions; if by hook or by crook we can spot two solutions, then we have 
the general solution. 

In particular, there is nothing to stop us trying to find a solution of the form 


for which 


x(t) = e rt 


• rt 

x = re 


(5.29) 


and 


*• 2 rt 

x = re 


Substituting into (5.28) we see that our guess (5.29) satisfies (5.28) if and only if 

r 2 + 2 fir + oo 2 = 0 (5.30) 


[an equation sometimes called the auxiliary equation for the differential equation 
(5.28)]. The solutions of this equation are, of course, r = ± y//3 2 — co\ Thus if 

we define the two constants 


r \ = ~ 




r 2 = 



(5.31) 


then the two functions e rit and e rit are two independent solutions of (5.28) and the 
general solution is 

x(t) = C x e r]t + C 2 e rit (5.32) 

= e~ pt (cyV/* 2 -"o' + o A . ( 5 . 33 ) 

This solution is rather too messy to be especially illuminating, but, by examining 
various ranges of the damping constant /?, we can begin to see what (5.33) entails. 


Undamped Oscillation 

If there is no damping then the damping constant f’< is zero, the square root in the 
exponents of (5.33) is just ico Q , and our solution reduces to 

x(t) = C x e ico + C 2 e~ ic ° ot , (5.34) 

the familiar solution for the undamped harmonic oscillator. 

Weak Damping 

Suppose next that the damping constant j J > is small. Specifically, suppose that 


P < co a , 


( 5 . 35 ) 
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a condition sometimes called underdamping. In this case, the square root in the 
exponents of (5.33) is again imaginary, and we can write 



where 

<jJ\ = yJ<o£ - P 2 . (5.36) 

The parameter co x is a frequency, which is less than the natural frequency co 0 . In the 
important case of very weak damping (/3 co 0 ), co x is very close to co Q . With this 

notation, the solution (5.33) becomes 

x(t) = e- p ' (cy'" 1 ' + C 2 e~ ilD A . (5.37) 

This solution is the product of two factors: The first, e ~^, is a decaying exponen¬ 
tial, which steadily decreases toward zero. The second factor has exactly the form 
(5.34) of undamped oscillations, except that the natural frequency co Q is replaced by 
the somewhat lower frequency co x . We can rewrite the second factor, as in Equation 
(5.11), in the form A cos (co x t — 8) and our solution becomes 

x(t) = Ae~^ cos(a) x t — 8). (5.38) 

This solution clearly describes simple harmonic motion of frequency co x with an 
exponentially decreasing amplitude Ae~^, as shown in Figure 5.11. The result (5.38) 
suggests another interpretation of the damping constant /3. Since j3 has the dimensions 
of inverse time, l/fi is a time, and we now see that it is the time in which the 
amplitude function Ae~^ falls to 1 /e of its initial value. Thus, at least for underdamped 
oscillations, can be seen as the decay parameter, a measure of the rate at which the 
motion dies out, 


(decay parameter) = fi [underdamped motion]. 



Figure 5.11 Underdamped oscillations can be thought of as sim¬ 
ple harmonic oscillations with an exponentially decreasing amplitude 
Ae~P f . The dashed curves are the envelopes, =t Ae~^ f . 
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The larger p the more rapidly the oscillations die out, at least for the case P < co 0 that 
we are discussing here. 


Strong Damping 

Suppose instead that the damping constant P is large. Specifically suppose that 

P > oo 0 , (5.39) 

a condition sometimes called overdamping. In this case, the square root in the 
exponents of (5.33) is real and our solution is 

x(t) = c+ C 2 « - ( /J+ ^ 2_a ’° 2 h (5.40) 

Here we have two real exponential functions, both of which decrease as time goes by 
(since the coefficients of t in both exponents are negative). In this case, the motion 
is so damped that it completes no bona fide oscillations. Figure 5.12 shows a typical 
case in which the oscillator was given a kick from O at t = 0; it slid out to a maximum 
displacement and then slid ever more slowly back again, returning to the origin only 
in the limit that t —> oo. The first term on the right of (5.40) decreases more slowly 
than the second, since the coefficient in its exponent is the smaller of the two. Thus 
the long-term motion is dominated by this first term. In particular, the rate at which 
the motion dies out can be characterized by the coefficient in the first exponent, 

(decay parameter) = ft — J ft 2 — co 2 [overdamped motion]. (5.41) 

Careful inspection of (5.41) shows that — contrary to what one might expect — the 
rate of decay of overdamped motion gets smaller if the damping constant P is made 
bigger. (See Problem 5.20.) 



Figure 5.12 Overdamped motion in which the oscillator is kicked 
from the origin at t = 0. It moves out to a maximum displacement and 
then moves back toward O asymptotically as t oo. 


Critical Damping 

The boundary between underdamping and overdamping is called critical damping 
and occurs when the damping constant is equal to the natural frequency, P = co 0 . This 
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case has some interesting features, especially from a mathematical point of view. 
When p = co 0 the two solutions that we found in (5.33) are the same solution, namely 

x(t) = e~ pt . (5.42) 

[This happened because the two solutions of the auxiliary equation (5.30) happen 
to coincide when ft — co 0 .\ This is the one case where our inspired guess, to seek a 
solution of the form x (/) = e rt , fails to find us two solutions of the equation of motion, 
and we have to find a second solution by some other method. Fortunately, in this case, 
it is not hard to spot a second solution: As you can easily check, the function 

x(t) = te~^ (5.43) 

is also a solution of the equation of motion (5.28) in the special case that ft = co 0 . (See 
Problems 5.21 and 5.24.) Thus the general solution for the case of critical damping is 

x{t) = C, e~ pt + C 2 1 e~ pt . (5.44) 

Notice that both terms contain the same exponential factor e~P‘. Since this factor is 
what dominates the decay of the oscillations as t -> oo, we can say that both terms 
decay at about the same rate, with decay parameter 

(decay parameter) = ft = co 0 [critical damping]. 

It is interesting to compare the rates at which the various types of damped oscil¬ 
lation die out. We have seen that in each case, this rate is determined by a “decay 
parameter,” which is just the coefficient of t in the exponent of the dominant expo¬ 
nential factor in x(t). Our findings can be summarized as follows: 


damping 

p 

decay parameter 

none 

o 

II 

0 

under 

P < <» 0 

p 

critical 

p=0 ),o 

p 

over 

P > 0) 0 

p-JP 2 -< 


Figure 5.13 is a plot of the decay parameter as a function of p and shows clearly that 
the motion dies out most quickly when p = co 0 ; that is, when the damping is critical. 
There are situations where one wants any oscillations to die out as quickly as possible. 
For example, one wants the needle of an analog meter (a voltmeter or pressure gauge, 
for instance) to settle down rapidly on the correct reading. Similarly, in a car, one 
wants the oscillations caused by a bumpy road to decay quickly. In such cases one 
must arrange for the oscillations to be damped (by the shock absorbers in a car), and 
for the quickest results the damping should be reasonably close to critical. 



Section 5.5 Driven Damped Oscillations 


179 



Figure 5.13 The decay parameter for damped oscillations as 
a function of the damping constant p. The decay parameter 
is biggest, and the motion dies out most quickly, for critical 
damping, with ft = co Q . 


5.5 Driven Damped Oscillations 


Any natural oscillator, left to itself, eventually comes to rest, as the inevitable damping 
forces drain its energy. Thus if one wants the oscillations to continue, one must arrange 
for some external “driving” force to maintain them. For example, the motion of the 
pendulum in a grandfather clock is driven by periodic pushes caused by the clock’s 
weights; the motion of a young child on a swing is maintained by periodic pushes 
from a parent. If we denote the external driving force by F(t ) and if we assume as 
before that the damping force has the form — bv, then the net force on the oscillator 
is —bv — kx + F(t) and the equation of motion can be written as 

mx + bx + kx = F(t). (5.45) 


Like its counterpart for undriven oscillations, this differential equation crops up in 
several other areas of physics. A prominent example is the LRC circuit of Figure 5.10. 
If we want the oscillating current in that circuit to persist, we must apply a driving 
EMF, £(0, in which case the equation of motion for the circuit becomes 


Lq + Rq + 


C 


q = £(/) 


(5.46) 


in perfect correspondence with (5.45). 

As before, we can tidy Equation (5.45) if we divide the equation by m and replace 
b/m by 2/3 and k/m by co ^. In addition, I shall denote F(t)/m by 


/« = 


F(t) 

m 


the force per unit mass. With this notation, (5.45) becomes 



(5.47) 
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Linear Differential Operators 

Before we discuss how to solve this equation, I would like to streamline our notation. 
It turns out that it is very helpful to think of the left side of (5.48) as the result of a 
certain operator acting on the function x(t). Specifically, we define the differential 
operator 

D =T-2 +1I5 T +CO o- (5-49) 

dt 2 dt 0 

The meaning of this definition is simply that when D acts on x it gives the left side 
of (5.48): 


Dx = x + 2 fix + CD^x. 

This definition is obviously a notational convenience — the equation (5.48) becomes 
just 


Dx = f (5.50) 

— but it is much more: The notion of an operator like (5.49) proves to be a powerful 
mathematical tool, with applications throughout physics. For the moment the impor¬ 
tant thing is that the operator is linear. We know from elementary calculus that the 
derivative of ax (where a is a constant) is just ax and that the derivative of x x + x 2 is 
just x x -f- x 2 . Since this also applies to second derivatives, it applies to the operator D : 

D(ax) = aDx and D(x x 4- x 2 ) = Dx x 4- Dx 2 . 

(Make sure you understand what these two equations mean.) We can combine these 
into a single equation: 


D(ax x 4- bx 2 ) = aDxi 4- bDx 2 (5.51) 

for any two constants a and b and any two functions Xj(t) and x 2 (t). Any operator 
that satisfies this equation is called a linear operator. 

We have actually used the property (5.51) of linear operators before. The equation 
(5.28) for a damped oscillator (not driven) can be written as 

Dx = 0. (5.52) 

The superposition principle asserts that if x x and x 2 are solutions of this equation, then 
so is ax x 4- bx 2 for any constants a and b. In our new operator notation, the proof is 
very simple: We are given that Dx x = 0 and Dx 2 = 0, and using (5.51) it immediately 
follows that 


D{ax x 4- bx 2 ) = aDx x 4- bDx 2 = 0 + 0 = 0; 
that is, ax x + bx 2 is also a solution. 

The equation (5.52), Dx — 0, for the undriven oscillator is called a homogeneous 
equation, since every term involves either x or one of its derivatives exactly once. The 
equation (5.50), Dx = /, is called an inhomogeneous equation, since it contains the 
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inhomogeneous term /, which does not involve v at all. Our job now is to solve this 
inhomogeneous equation. 


Particular and Homogeneous Solutions 

Using our new operator notation, we can find the general solution of Equation (5.48) 
surprisingly easily; in fact, we’ve already done most of the work. Suppose first that we 
have somehow spotted a solution; that is, we’ve found a function x p (t) that satisfies 

Dx p - /. (5.53) 

We call this function x p (t) a particular solution of the equation, and the subscript “p” 
stands for “particular.” Next let us consider for a moment the homogeneous equation 
Dx = 0 and suppose we have a solution jt h (r), satisfying 

Dx h = 0. (5.54) 

We’ll call this function a homogeneous solution. 8 We already know all about the 
solutions of the homogeneous equation, and we know from (5.32) that x h must have 
the form 


jc h (0 = C x e rit + <V 2 ', (5.55) 

where both exponentials die out as t -> oo. 

We’re now ready to prove the crucial result. First, if x p is a particular solution 
satisfying (5.53), then x p + x h is another solution, for 

D(x p + x^) = Dx p + Dxh — f + 0 = /. 

Given the one particular solution x p , this gives us a large number of other solutions 
x p + x h . And we have in fact found all the solutions, since the function x h contains 
two arbitrary constants, and we know that the general solution of any second-order 
equation contains exactly two arbitrary constants. Therefore x p + jc h , with x h given 
by (5.55) is the general solution. 

This result means that all we have to do is somehow to find a single particular 
solution x p (t) of the equation of motion (5.48), and we shall have every solution in 
the form x(t) = x p (t) + x h (t). 


Complex Solutions for a Sinusoidal Driving Force 

I shall now specialize to the case that the driving force f(t) is a sinusoidal function 
of time, 


fit) = f 0 cos(cut). 


(5.56) 


8 Another common name is the complementary function. This has the disadvantage that it’s hard 
to remember which is “particular” and which “complementary ” 
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where f 0 denotes the amplitude of the driving force [actually the amplitude divided 
by the oscillator’s mass, since f(t) = F(t)/m] and co is the angular frequency of 
the driving force. (Be careful to distinguish the driving frequency co from the natural 
frequency co Q of the oscillator. These are entirely independent frequencies, although 
we shall see that the oscillator responds most when co ^ co 0 .) The driving force for 
many driven oscillators is at least approximately sinusoidal. For example, even the 
parent pushing the child on a swing can be crudely approximated by (5.56); the driving 
EMF induced in your radio’s circuits by a broadcast signal is almost perfectly of this 
form. Probably the chief importance of sinusoidal driving forces is that, according 
to Fourier’s theorem, 9 essentially any driving force can be built up as a series of 
sinusoidal forces. 

Let us therefore assume the driving force is given by (5.56), so that the equation 
of motion (5.48) takes the form 

x + 2fx + co^x = f Q cos(cot). (5.57) 

Solving this equation is greatly simplified by the following trick: For any solution of 
(5.57), there must be a solution of the same equation, but with the cosine on the right 
replaced by a sine function. (After all, these two differ only by a shift in the origin of 
time.) Accordingly, there must also be a function y(t ) that satisfies 

y + 2 fiy + co 2 y = f 0 sin(tu?). (5.58) 

Suppose now we define the complex function 

z(t) =x(t) 4- iy(t), (5.59) 

with x(t) as its real part and y (t) as its imaginary part. If we multiply (5.58) by i and 
add it to (5.57), we find that 


z + 2 f>z + co^z — f 0 e lcot . (5.60) 

Although it may not yet look it, Equation (5.60) is a tremendous advance. Because 
of the simple properties of the exponential function, (5.60) is much easier to solve 
than either (5.57) or (5.58). And as soon as we find a solution z(t) of (5.60), we have 
only to take its real part to have a solution of the equation (5.57) whose solutions we 
actually want. 

In seeking a solution of (5.60), we are obviously free to try any function we like. 
In particular, let’s see if there is a solution of the form 

z(t) = Ce iQ)t , (5.61) 

where C is an as yet undetermined constant. If we substitute this guess into the left 
side of (5.60), we get 


(-co 2 + 2 ipco + oj 2 )Ce imt = f 0 e iwt . 


9 Named for the French mathematician, the Baron Jean Baptiste Joseph Fourier, 1768-1830. See 
Sections 5.7-5.9. 
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In other words, our guess (5.61) is a solution of (5.60) if and only if 


C = 


fo 

co 2 — co 2 + 2i/3co ’ 


(5.62) 


and we have succeeded in finding a particular solution of the equation of motion. 

Before we take the real part of z(t) = Ce ia)t , it is convenient to write the complex 
coefficient C in the form 

C = Ae~ iS , (5.63) 


where A and 8 are real. [Any complex number can be written in this form; the particular 
notation is chosen to match (5.14).] To identify A and 8 , we must compare (5.62) and 
(5.63). First, taking the absolute value squared of both equations we find that 


A 2 = CC* = 


fo 

co 2 — co 2 + 2 i/3co 


fo 

0)2 — co 2 — 2i/3co 


or 




(co 2 — co 2 ) 2 + 4 f3 2 co 2 


(5.64) 


(Make sure you understand this derivation. See Problem 5.35 for some guidance.) 
We are going to see in a moment that A is just the amplitude of the oscillations 
caused by the driving force f(t). Thus the result (5.64) is the most important result 
of this discussion. It shows how the amplitude of oscillations depends on the various 
parameters. In particular, we see that the amplitude is biggest when co 0 & co, so that 
the denominator is small; in other words, the oscillator responds best when driven at 
a frequency co that is close to its natural frequency co Q , as you would probably have 
guessed. 

Before we continue to discuss the properties of our solution, we need to identify 
the phase angle 8 in (5.63). Comparing (5.63) and (5.62) and rearranging, we see that 


f Q e l8 — A(co 2 — co 2 + 2i/3oo). 


Since f Q and A are real, this says that the phase angle 8 is the same as the phase 
angle of the complex number (co 2 — co 2 ) + 2 ifioo. This relationship is illustrated in 
Figure 5.14, from which we conclude that 


8 = arctan 


2/3 co 
co 2 — co 2 

o 



(5.65) 


Our quest for a particular solution is now complete. The “fictitious” complex 
solution introduced in (5.59) is 


z(t) = Ce ia)t = Ae i(Mt ~ 8) 


and the real part of this is the solution we are seeking, 

x(t) = Acos(cot — 8) (5.66) 

where the real constants A and 8 are given by (5.64) and (5.65). 



184 


Chapter 5 Oscillations 



Figure 5.14 The phase angle 8 is the angle of this triangle. 


The solution (5.66) is just one particular solution of the equation of motion. The 
general solution is found by adding any solution of the corresponding homogeneous 
equation, as given by (5.55); that is, the general solution is 

x(t) = A cos (cot - 8) + C x e rit + C 2 e rit . (5.67) 

Because the two extra terms in this general solution both die out exponentially as time 
passes, they are called transients. They depend on the initial conditions of the problem 
but are eventually irrelevant: The long-term behavior of our solution is dominated by 
the cosine term. Thus the particular solution (5.66) is the solution with which we are 
usually concerned, and we shall explore its properties in the next section. 

Before we discuss an example of the motion (5.67), it is important that you be 
very clear as to the type of system to which (5.67) applies, namely, any oscillator 
for which both the restoring force (—kx) and the resistive force (~bx) are linear — 
a driven, damped linear oscillator, whose equation of motion (5.45) is a linear 
differential equation. Because nonlinear differential equations are often hard to solve, 
most mechanics texts have until recently focussed on linear equations. This created 
the false impression that linear equations were in some sense the norm, and that the 
solution (5.67) was the only (or, at least, the only important) way for an oscillator to 
behave. As we shall see in Chapter 12 on nonlinear mechanics and chaos, an oscillator 
whose equation of motion is nonlinear can behave in ways that are astonishingly 
different from (5.67). One important reason for studying the linear oscillator here 
is to give you a backdrop against which to study the nonlinear oscillator later on. 

The details of the motion (5.67) depend on the strength of the damping param¬ 
eter To be specific, let us assume that our oscillator is weakly damped, with less 
than the natural frequency co 0 (underdamping). In this case, we know that the two 
transient terms of (5.67) can be rewritten as in (5.38), so that 



You need to think very carefully about this potentially confusing formula. The second 
term on the right is the homogeneous or transient term, and I have added the subscript 
“tr” to distinguish the constants A tr and <$ tr from the A and 8 of the first term. The 
two constants A tr and <5 tr are arbitrary constants; (5.68) is a possible motion of our 
system for any values of A tr and <$ tr , which are determined by the initial conditions. 
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The factor makes clear that this transient term decays exponentially and is indeed 
irrelevant to the long-term behavior. As it decays, the transient term oscillates with 
the angular frequency co l of the undriven (but still damped) oscillator, as in (5.36). 
The first term is our particular solution and its two constants A and 8 are certainly 
not arbitrary; they are determined by (5.64) and (5.65) in terms of the parameters of 
the system. This term oscillates with the frequency co of the driving force and with 
unchanging amplitude, for as long as the driving force is maintained. 

example 5.3 Graphing a Driven Damped Linear Oscillator 

Make a plot of x (t) as given by (5.68) for a driven damped linear oscillator which j 
is released from rest at the origin at time t = 0, with the following parameters: 
Drive freqency co = 2n, natural frequency co 0 = 5co, decay constant fi = co 0 / 20, 
and driving amplitude f Q = 1000. Show the first five drive cycles. 

The choice of drive frequency equal to 2ic means that the drive period is J 
r = 2tc I co — 1; this means simply that we have chosen to measure time in units | 
of the drive period — often a convenient choice. That co 0 = 5 co means that our 
oscillator has a natural frequency five times the drive frequency; this will let j 
us distinguish easily between the two on a graph. That fi = coj 20 means that j 
the oscillator is rather weakly damped. Finally the choice of f 0 — 1000 is just j 
a choice of our unit of force; the reason for this odd-seeming choice is that it | 
leads to a conveniently sized amplitude of oscillation (namely A close to 1). 

Our first task is to determine the various constants in (5.68) in terms of the 
given parameters. In fact this is easier if we rewrite the transient term of (5.68) j 
in the “cosine plus sine” form, so that 

x(t) — Acos(cot — 8) + e~~^[BiCos(coit) + B 2 sin(co l t)]. (5.69) | 

| 

j 

The constants A and 8 are determined by (5.64) and (5.65), which, for the given 
parmeters, yield 


A = 1.06 and 8 = 0.0208. 

The frequency co 1 is 

co, = Ja>l - p 2 = 9Mln, 

which is very close to co 0 , as we would expect for a weakly damped oscillator. 
To find B x and B 2 , we must equate x(0) as given by (5.69) to its given initial 
value x 0 , and likewise the corresponding expression for x 0 to the initial value v 0 . 
This gives two simultaneous equations for B x and B 2 , which are easily solved 
(Problem 5.33) to give 

1 

B\ = x 0 ~ A cos 8 and B 2 — — (y 0 — coA sin <5 + /3B{) (5.70) 

<x>i 

or, with the numbers, including the initial conditions x 0 = v 0 = 0, 

B x = —1.05 and B 2 = -0.0572. 
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Figure 5.15 The response of a damped, linear oscillator to a si¬ 
nusoidal driving force, with the time t shown in units of the drive 
period, (a) The driving force is a pure cosine as a function of time, 
(b) The resulting motion for the initial conditons x 0 = v Q = 0. For 
the first two or three drive cycles, the transient motion is clearly vis¬ 
ible, but after that only the long-term motion remains, oscillating 
sinusoidally at exactly the drive frequency. As explained in the text, 
the sinusoidal motion after t & 3 is called an attractor. 


\ 


Putting all of these numbers into (5.69), we can now plot the motion, as 

I 

shown in Figure 5.15, where part (a) shows the driving force f(t ) = f 0 cos (cot) 

| and part (b) the resulting motion ;t(0 of the oscillator. The driving force is, of 

| 

course, perfectly sinusoidal with period 1. The resulting motion is much more 
interesting. After about three drive cycles (t ^ 3), the motion is indistinguishable 
from a pure cosine, oscillating at exactly the drive frequency; that is, the 
transients have died out and only the long-term motion remains. Before t ^ 3, 

I 

however, the effects of the transients are clearly visible. Since they oscillate at 
the faster natural frequency co 0 , they show up as a rapid succession of bumps 
and dips. In fact, you can easily see that there are five such bumps within the 
first drive cycle, indicating that co 0 = 5co. 


Because the transient motion depends on the initial values x 0 and v 0 , different 
values of x 0 and v 0 would lead to quite different initial motion. (See Problem 5.36.) 
After a short time however (a couple of cycles in this example), the initial differences 
disappear and the motion settles down to the same sinusoidal motion of the particular 
solution (5.66), irrespective of the initial conditions. For this reason, the motion 
of (5.66) is sometimes called an attractor — the motions corresponding to several 
different initial conditions are “attracted” to the particular motion (5.66). For the linear 
oscillator discussed here, there is a unique attractor (for a given driving force): Every 
possible motion of the system, whatever its initial conditions, is attracted to the same 
motion (5.66). We shall see in Chapter 12 that for nonlinear oscillators there can be 
several different attractors and that for some values of the parameters the motion of 
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an attractor can be far more complicated than simple harmonic oscillation at the drive 
frequency. 

The amplitude and phase of the attractor seen in Figure 5.15(b) depend on the 
parameters of the driving force (but not, of course, on the initial conditions). The 
dependence of the amplitude and phase on these parameters is the subject of the next 
section. 


5.6 Resonance 


In the previous section we considered a damped oscillator that is being driven by 
a sinusoidal driving force (actually force divided by mass) fit) = f Q cos(cot ) with 
angular frequency co. We saw that, apart from transient motions that die out quickly, 
the system’s response is to oscillate sinusoidally at the same frequency, co: 

x(t) — A cos (cot — 8 ), 


with amplitude A given by (5.64), 


A 


2 



f 2 

Jo _ 

co 2 ) 2 + 4/3 2 co 2 ’ 


(5.71) 


and phase shift 8 given by (5.65). 

The most obvious feature of (5.71) is that the amplitude A of the response is 
proportional to the amplitude of the driving force, A oc / G , a result you might have 
guessed. More interesting is the dependence of A on the frequencies co Q (the natural 
frequency of the oscillator) and co (the frequency of the driver), and on the damping 
constant /?. The most interesting case is when the damping constant f> is very small, 
and this is the case I shall discuss. With /3 small, the second term in the denominator of 
(5.71) is small. If co 0 and co are very different, then the first term in the denominator of 
(5.71) is large, and the amplitude of the driven oscillations is small. On the other hand, 
if co 0 is very close to co , both terms in the denominator are small, and the amplitude 
A is large. This means that if we vary either co 0 or co , there can be quite dramatic 
changes in the amplitude of the oscillator’s motion. This is illustrated in Figure 5.16, 
which shows A 2 as a function of co 0 with co fixed, for a rather weakly damped system 
i/3 = O.ko). (Note that, because the energy of the system is proportional to A 2 , it is 
usual to make plots of A 2 rather than A.) 

Although the behavior illustrated in Figure 5.16 is startlingly dramatic, the qualita¬ 
tive features are what you might have expected. Left to its own devices, the oscillator 
vibrates at its natural frequency co 0 (or at the slightly lower frequency co 1 if we allow 
for the damping). If we try to force it to vibrate at a frequency co , then for values of 
co close to co 0 the oscillator responds very well, but if co is far from co 0 , it hardly re¬ 
sponds at all. We refer to this phenomenon — the dramatically greater response of an 
oscillator when driven at the right frequency — as resonance. 

An everyday application of resonance is the reception of radio signals by the LRC 
circuit in your radio. As we saw, the equation of motion of an LRC circuit is exactly 
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Figure 5.16 The amplitude squared, A 2 , of a driven oscillator, 
shown as a function of the natural frequency co 0 , with the driving 
frequency co fixed. The response is dramatically largest when co Q 
and co are close. 


the same as that of a driven oscillator, and LRC circuits show the same phenomenon 
of resonance. When you tune your radio to receive station KVOD at 90.1 MHz, you 
are adjusting an LRC circuit in the radio so that its natural frequency is 90.1 MHz. 
The many radio stations in your neighborhood are all sending out signals, each at its 
own frequency and each inducing a tiny EMF in the circuit of your radio. But only 
the signal with the right frequency actually succeeds in driving an appreciable current, 
which mimics the signal sent out by your favorite KVOD and reproduces its broadcast 
sounds. 

An example of a mechanical resonance of the kind discussed here is the behavior 
of a car driving on a “washboard” road that has worn into a series of regularly spaced 
bumps. Each time a wheel crosses a bump, it is given an upward impulse, and the 
frequency of these impulses depends on the car’s speed. There is one speed at which 
the frequency of these impulses equals the natural frequency of the wheels’ vibration 
on the springs 10 and the wheels resonate, causing an uncomfortable ride. If the car 
drives slower or faster than this speed, it goes “off resonance” and the ride is much 
smoother. 

Another example occurs when a platoon of soldiers marches across a bridge. A 
bridge, like almost any mechanical system, has certain natural frequencies, and if the 
soldiers happened to march with a frequency equal to one of these natural frequencies, 
the bridge could conceivably resonate sufficiently violently to break. For this reason, 
soldiers break step when marching across a bridge. 

The details of the resonance phenomenon are a bit complicated. For example, the 
exact location of the maximum response depends on whether we vary co 0 with co fixed 
or vice versa. The amplitude A is a maximum when the denominator, 

denominator = ( co 2 — co 2 ) 2 + Af5 2 co 2 (5.72) 


10 It is the wheels (plus axles) that exhibit the resonant oscillations; the much heavier body of 
the car is relatively unaffected. 
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of (5.71) is a minimum. If we vary co 0 with co fixed (as when you tune a radio to pick 
up your favorite station) this minimum obviously occurs when co 0 = co , making the 
first term zero. On the other hand, if we vary co with co Q fixed (which is what happens 
in many applications), then the second term in (5.72) also varies and a straightforward 
differentiation shows that the maximum occurs when 

60 = co 2 = J 0*1 - 2£ 2 . (5.73) 

However, when /3 <$C co 0 (as is usually the most interesting case), the difference 
between (5.73) and co = co 0 is negligible. 

We have met so many different frequencies in this chapter that it may be worth 
pausing to review them. First there is co 0 the natural frequency of the oscillator (in 
the absence of any damping). Next when we added in a little damping, we found 

that the same system oscillated sinusoidally with frequency co l = y4o 2 — /3 2 under an 
exponentially decaying envelope. Then we added a driving force with frequency co, 
which can, in principle, take on any value independently of the previous two. However, 
the response of the driven oscillator is biggest when co & co 0 ; specifically, if we vary 
co with co Q fixed, the maximum response comes when co = co 2 as defined by (5.73). To 
summarize: 

co 0 = y/k/m — natural frequency of undamped oscillator 

co i — _ p 2 __ frequency of damped oscillator 

co — frequency of driving force 

co 2 = tJco^ — 2fi 2 = value of co at which response is maximum. 

In any case, the maximum amplitude of the driven oscillations is found by putting 
co 0 & co in (5.71), to give 


A 


max 


fo 

2 fico 0 


(5.74) 


This shows that smaller values of the damping constant /3 lead to larger values of the 
maximum amplitude of oscillation, as illustrated in Figure 5.17. 11 


Width of the Resonance; the Q Factor 

You can see clearly from Figure 5.17 that if we make the damping constant f3 smaller, 
the resonance peak not only gets higher, but also gets narrower. We can make this idea 
more precise by defining the width (or full width at half maximum or FWHM) as 
the interval between the two points where A 2 is equal to half its maximum height. It 


11 In this figure I chose to plot A 2 against co with co 0 fixed, rather than the other way around as 
in Figure 5.16. Note that the curves have very similar shapes either way. 
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Figure 5.17 The amplitude for driven oscillations as a function 
of the driving frequency co for three different values of the damp¬ 
ing constant p. Note how as ft decreases the resonance peak gets 
higher and sharper. 


is a simple exercise (Problem 5.41) to show that the two half-maximum points are at 
co & co 0 ± /3, as in Figure 5.18. Thus the full width at half maximum is 

FWHM ^ 2/3 (5.75) 


or, equivalently, the half width at half maximum is 

HWHM « p. (5.76) 


The sharpness of the resonance peak is indicated by the ratio of its width 2ft to its 
position, co Q . For many purposes, we want a very sharp resonance, so it is common 
practice to define a quality factor Q as the reciprocal of this ratio: 


Q = 



(5.77) 



Figure 5.18 The full width at half maximum (FWHM) is the 
distance between the points where A 2 is half its maximum value. 
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A large Q indicates a narrow resonance, and vice versa. For example, clocks depend 
on the resonance in an oscillator (a pendulum, or a quartz crystal, for instance) to 
regulate the mechanism to move with very well-defined frequency. This requires that 
the width 2/3 be very small compared to the natural frequency co Q . In other words, a 
good clock needs a high Q. The Q for a typical pendulum may be around 100; that 
for a quartz crystal around 10,000. Therefore quartz watches keep much better time 
than a typical grandfather clock. 12 

There is another way to look at the quality factor Q. We saw that, in the absence 
of any driving force, the oscillations die out in a time of order 1//3, 


(decay time) = 1//3. 

(This was actually the time for the amplitude to drop to l/e of its initial value.) The 
period of a single oscillation is, of course, 


period = 2?r/co Q . 


(Remember we’re assuming that co 0 , so we don’t need to distinguish between co 0 
and co Thus we can rewrite the definition of Q as 


co Q l/p decay time 

Q = — = 7T -= 7r- 

2fi 2 tc/co 0 period 


(5.78) 


The ratio on the right is just the number of periods in the decay time. Thus, the quality 
factor Q is Jt times the number of cycles our oscillator makes in one decay time. 13 


The Phase at Resonance 

The phase difference 8 by which the oscillator’s motion lags behind the driving force 
is given by (5.65) as 


8 = arctan ( ) . (5.79) 

W-^7 

Let us follow this phase as we vary co , starting well below a narrow resonance (/3 
small). With co <$C &> 0 , (5.79) implies that <5 is very small; that is, while co <^co 0 the 
oscillations are almost perfectly in step with the driving force. (This was the case 
in Figure 5.15.) As co is increased toward co 0 , so 8 slowly increases. At resonance, 
where co = co 0 , the argument of the arctangent in (5.79) is infinite, so 8 = it/2 and 
the oscillations are 90° behind the driving force. Once co > co 0 , the argument of the 


12 Actually, both quartz watches and grandfather clocks keep much better time than this simple 
discussion would suggest. A good chronometer keeps the frequency very close to the center of 
the resonance. Thus the variability of the frequency is actually much smaller than the width of the 
resonance. Nevertheless, the stated conclusion is correct. 

13 Yet another definition (and perhaps the most fundamental) is that Q = 2n times the ratio of 
the energy stored in the oscillator to the energy dissipated in one cycle. See Problem 5.44. 
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Figure 5.19 The phase shift 8 increases from 0 through 7r/2 
to 7i as the driving frequency co passes through resonance. The 
narrower the resonance, the more suddenly this increase occurs. 
The solid curve is for a relatively narrow resonance (/? = 0.03&> o 
or Q — 16.7), and the dashed curve is for a wider resonance 
00 = O3co 0 or Q = 1.67). 


arctangent is negative and approaches 0 as co increases; thus 8 increases beyond n/2 
and eventually approaches jr. In particular, once co co 0 , the oscillations are almost 
perfectly out of step with the driving force. All of this behavior is illustrated for 
two different values of f3 in Figure 5.19. Notice, in particular, that the narrower the 
resonance, the more quickly 8 jumps from 0 to tt. 

In the resonances of classical mechanics, the behavior of the phase (as in Fig¬ 
ure 5.19) is usually less important than that of the amplitude (as in Figure 5.18). 14 In 
atomic and nuclear collisions, the phase shift is often the quantity of primary inter¬ 
est. Such collisions are governed by quantum mechanics, but there is a corresponding 
phenomenon of resonance. A beam of neutrons, for example, can “drive” a target 
nucleus. When the energy of the beam equals a resonant energy of the system (in 
quantum mechanics energy plays the role of frequency) a resonance occurs and the 
phase shift increases rapidly from 0 to tt. 


5.7 Fourier Series* 


* Fourier series have broad application in almost every area of modern physics. Nevertheless, 
we shall not be using them again until Chapter 16. Thus, if you are pressed for time you could 
omit the last three sections of this chapter on a first reading. 

In the last two sections, we have discussed an oscillator that is driven by a sinusoidal 
driving force f(t) = f 0 cos(cot). There are two main reasons for the importance of 


14 The behavior of 8 can, nevertheless, be observed. Make a simple pendulum from a piece of 
string and a metal nut, and drive it by holding it at the top and moving your hand from side to side. 
The most obvious thing is that you will be most successful at driving it when your frequency equals 
the natural frequency, but you can also see that when you drive more slowly the pendulum moves 
in step with your hand, whereas when you move more quickly the pendulum moves oppositely to 
your hand. 
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Figure 5.20 Two examples of periodic functions with period r. (a) A rectangular 
pulse, which could represent a hammer hitting a nail with a constant force at 
intervals of r, or a digital signal in a telephone line, (b) A smooth periodic signal, 
which could be the pressure variation of a musical instrument. 


sinusoidal driving forces: The first is simply that there are many important systems 
in which the driving force is sinusoidal — the electrical circuit in a radio is a good 
example. The second is somewhat subtler. It turns out that any periodic driving force 
can be built up from sinusoidal forces using the powerful technique of Fourier series. 
Thus, in a sense that I shall try to describe, by solving the motion with a sinusoidal 
driver we have already solved the motion with any periodic driver. Before we can 
appreciate this wonderful result, we need to review some aspects of Fourier series. In 
this section I sketch the needed properties of Fourier series; 15 in the next we can apply 
them to the driven oscillator. 

Let us consider a function f(t) that is periodic with period r; that is, the function 
repeats itself every time t advances by the period r: 

fit + r) = fit) 

whatever the value of t. We can describe a function with this property as being r- 
periodic. A simple example of a r-periodic function would be the force exerted on a 
nail by a hammer that is being swung at intervals of r, as sketched in Figure 5.20(a). 
Another could be the pressure exerted on your ear drum by a note played by a musical 
instrument, as sketched in 5.20(b). It is easy to think up many more periodic functions. 
In particular, there are lots of sinusoidal functions that are periodic with any given 
period: The functions 

cosilnt/x), cos(47zT/r), cos(6jrr/r), ••• (5.80) 

are all r-periodic, as are the corresponding sine functions. (If t is increased by the 
amount r, each of these functions returns to its original value — see Figure 5.21.) 
We can write these sinusoidal functions a little more compactly if we introduce the 
angular frequency co = 2n/r, in which case all of the functions of (5.80) and the 
corresponding sines can be written as 

cos (ncot) and sin (ncot) [n = 0, 1, 2, • • •]. (5.81) 


15 As usual, I shall try to describe all of the theory that we shall be needing. For more details, see, 
for example, Mathematical Methods in the Physical Sciences by Mary Boas (Wiley, 1983), Ch. 7. 
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Figure 5.21 Any function of the form cos(2njtt/z) (or the corresponding sine) 
is periodic with period r if n is an integer. Notice that cos(47rt/r) also has the 
smaller period r/2, but this doesn’t change the fact that it has period r as well. 


(If n = 0 the cosine function is just the constant 1 — which is certainly periodic — 
while the sine is 0 and not at all interesting.) 

That the sine and cosine functions (5.81) are all r-periodic is reasonably obvious. 
(Be sure you can see this.) What is truly amazing is that, in a sense, these sine and 
cosine functions define all possible r-periodic functions: In 1807 the French math¬ 
ematician Jean Baptiste Fourier (1768-1830) realized that every r-periodic function 
can be written as a linear combination of the sines and cosines of (5.81). That is, if 
fit) is any 16 periodic function with period r then it can be expressed as the sum 


fit) = yy la„cos(nm) + h„ sin(»a>0] 




where the constants a n and b n depend on the function fit). This extraordinarily useful 
result is called Fourier’s theorem, and the sum (5.82) is called the Fourier series for 
fit). 

It is not hard to see why Fourier’s theorem met with considerable surprise, and 
even skepticism, when he first published it. It claims that a discontinous function, 
such as the rectangular pulse of Figure 5.20(a), can be built up with sine and cosine 
functions that are continuous and perfectly smooth. Surprising or not, this turns out 
to be true, as we shall see by example shortly. Perhaps even more surprising, it is 
often the case that one gets an excellent approximation by retaining just the first few 
terms of a Fourier series. Thus, instead of having to handle a fairly disagreeable and 
possibly discontinuous function, we have only to handle a reasonably small number 
of sines and cosines. Before we discuss the application of Fourier’s theorem to the 
driven oscillator, we need to look at a few properties of Fourier series. 

The proof of Fourier’s theorem is difficult — indeed it was many years after 
Fourier’s discovery before a satisfactory proof was found — and I shall simply ask 
you to accept it. However, once the result is accepted it is easy to learn to use it. In 


16 As always with theorems of this kind, there are certain restrictions on the “reasonableness” 
of the function fit), but certainly Fourier’s theorem is valid for all of the functions we shall have 
occasion to use. 
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particular, for any given periodic function f(t) it is easy to find the coefficients a n 
and b n . Problem 5.48 gives you the opportunity to show that these coefficients are 
given by 



and 



Unfortunately the coefficients for n = 0 require separate attention. Since the term 
sin ncot in (5.82) is identically zero for n — 0, the coefficient b 0 is irrelevant and we 
can simply define it to be zero. It is very easy to show (Problem 5.46) that 



Armed with these formulas for the Fourier coefficients, it is easy to find the Fourier 
series for any given periodic function. In the following example, we do this for the 
rectangular pulses of Figure 5.20(a). 




example 5.4 Fourier Series for the Rectangular Pulse 


Find the Fourier series for the periodic rectangular pulse /(f) shown in Figure 
5.22 in terms of the period t, the pulse height / max , and the duration of the pulse 
At. Using the values r = 1, / max = 1, and At = 0.25, plot /(f), as well as the 
sum of the first three terms of its Fourier series, and the sum of the first eleven 
terms. 

Our first task is to calculate the Fourier coefficients a n and b n for the given 
function. First according to (5.85) the constant term a 0 is 


1 

a 0 = - 

T 

_ 1 
X 



f At 

J max L 


X 


(5.86) 
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Figure 5.22 A periodic rectangular pulse. The period is r, the 
duration of the pulse is At, and the pulse height is / max . 


where the change in the limits of integration was allowed because the integrand 
fit) is zero outside ± Ar/2. Next, according to (5.83), all the others coefficients 
(n > 1) are given by 


I 


2 f r/2 

a n = — / f(t)cos(ncot)dt 

r J—z/2 


2 /, 


max 


/ 


At12 


cos {ncot) dt 


—Ar/2 


4/ max r Ar/2 (2nnt\ . 2/ max . (nn Ax\ n _. 

= __7max / cos - ^ Sin - . (5.87) 

r Jo \ r ) rcn \ x ) 

I 

| 

Notice that in passing from the second to the third line, I used a trick that is 
often useful in evaluating Fourier coefficients. The integrand on the second line, 
cos (ncot), is an even function; that is, it has the same value at any point t as at 
—t. Therefore we could replace any integral from — T to T by twice the integral 
from 0 to T. 

Finally the b coefficients are all exactly zero, for if you examine the integral 
(5.84) you will see that (in this case) the integrand is an odd function; that is, its 
value at any point t is the negative of its value at —t. [Moving from t to—t leaves 
f(t ) unchanged, but reverses the sign of sin (ncot).] Therefore any integral from 
— T to T is zero, since the left half exactly cancels the right. 

The required Fourier series is therefore 


| 00 

j f(t) = #0 + E a n cosincot) (5.88) 

| n =1 

with the constant term a 0 given by (5.86) and all the remaining a coefficients 
{n > 1) by (5.87). If we put in the given numbers, these coefficients can all be 

I 

evaluated and the resulting Fourier series is 

i ^ 

i 

fit) = f m ax [0.25 + 0.24 cos(27tO + 0.23cos(47rt) + 0.20cos(67rt) 

I 

| + 0.16cos(8jrt) + 0.12cos(107rt) + 0.08cos(127rf) + • • • ] (5.89) 
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Figure 5.23 (a) The sum of the first three terms of the Fourier series for the 

rectangular pulse of Figure 5.22. (b) The sum of the first 11 terms. 


The practical value of a Fourier series is usually greatest if the series con¬ 
verges rapidly, so that we get a reliable approximation by retaining just the first 
few terms of the series. Figure 5.23(a) shows the sum of the first three terms of 
the series (5.89) and the rectangular pulse itself. As you might expect, with just 
three smooth terms we do not get a sensationally accurate approximation to the 
original discontinuous function. Nevertheless, the three terms do a remarkable 
job of imitating the general shape. By the time we have included 11 terms, as in 
Figure 5.23(b), the fit is amazing. 17 In the next section, we shall use the method 
of Fourier series to solve for the motion of an oscillator driven by the periodic 
pulses of this example. We shall find the solution as a Fourier series, which con¬ 
verges so quickly that just the first 3 or 4 terms tell us most of what is interesting 
to know. 


5.8 Fourier Series Solution for the Driven Oscillator* 


* This section contains a beautiful application of the method of Fourier series. Important as it 
is to understand this method, you can nevertheless omit the section without loss of continuity. 

In this section we shall combine our knowledge of Fourier series (Section 5.7) with 
our solution of the sinusoidally driven oscillator (Section 5.5) to solve for the motion 
of an oscillator that is driven by an arbitrary periodic driving force. To see how this 
works, let us return to the equation of motion (5.48) 

x + 2 fix + oo^x = / 


where x = x{t) is the position of the oscillator, f3 is the damping constant, co 0 is the 
natural frequency, and f — f{t) is any periodic driving force (actually force/mass) 
with period r. As before it is convenient to rewrite this in the compact form 

Dx = f 


17 Notice, however, that the Fourier series still has a little difficulty in the immediate neigh¬ 
borhood of the discontinuities in f(t). This tendency for the Fourier series to overshoot near a 
discontinuity is called the Gibbs phenomenon. 
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where D stands for the linear differential operator 


D = + 2j6— + co 2 

dt 2 dt 0 


The use of Fourier series to solve this problem hinges on the following observation: 
Suppose that the force fit) is the sum of two forces, fit) = f x (t) + f 2 it), for each 
of which we have already solved the equation of motion. That is, we already know 
functions x x (t) and x 2 (t) that satisfy 


Dx x = /j and Dx 2 = / 2 . 

Then the solution 18 to the problem of interest is just the sumi(f) = x x (t) + x 2 (t ), as 
we can easily show: 


Dx = D(x { + x 2 ) = Dx 1 + Dx 2 = /i + f 2 = f 

where the crucial second step is valid because D is linear. This argument would work 
equally well however many terms were in the sum for fit), so we have the conclusion: 
If the driving force fit) can be expressed as the sum of any number of terms 

/(0 = £/,(*) 

n 

and if we know the solutions x n (t ) for each of the individual forces f n (t) 9 then the 
solution for the total driving force fit) is just the sum 

x(t) = Y^x n (t). 

n 

This result is ideally suited for use in combination with Fourier’s theorem. Any 
periodic driving force fit) can be expanded in a Fourier series of sines and cosines, 
and we already know the solutions for sinusoidal driving forces. Thus, by adding these 
sinusoidal solutions together, we can find the solution for any periodic driving force. 
To simplify our writing, let us suppose that the driving force f(t ) contains only cosine 
terms in its Fourier series. [This was the case for the rectangular pulse of Example 5.4, 
and is true for any even function — satisfying f(—t ) = fit) — because this condition 
guarantees that the coefficients of the sine terms are all zero.] In this case, the driving 
force can be written as 



18 Strictly speaking we should not speak of the solution, since our second-order differential 
equation has many solutions. However, we know that the difference between any two solutions is 
transient — decays to zero — and our main interest is in the long-term behavior, which is therefore 
essentially unique. 
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where f n denotes the nth Fourier coefficient of /(f), and co = lit/x as usual. Now, 
each individual term f n cos {ncot) has the same form (5.56) that we assumed for the 
sinusoidal driving force in Section 5.5 (except that the amplitude f 0 has become f n and 
the frequency co has become nco). The corresponding solution was given in (5.66), 19 

x n (t) = A n cos (ncot — 8 n ) (5.91) 


where 



_ fn _ 

J (&>% — n 2 o> 2 ) 2 + 4/3 2 n 2 <w 2 


from (5.64), and 


(5.92) 


8 n = arctan 



2 finco 


2 — n 2 
o 


CO 



(5.93) 


from (5.65). Since (5.91) is the solution for the driving force f n cos (ncot), the solution 
for the complete force (5.90) is the sum 


00 

x (f) = A n cos (ncot — 8 n ). (5,94) 

0 


This completes the solution for the long-term motion of an oscillator driven by a 
periodic driving force /(f). To summarize, the steps are: 

1. Find the coefficients f n in the Fourier series (5.90) for the given driving force 
fit). 

2. Calculate the quantities A n and 8 n as given by (5.92) and (5.93). 

3. Write down the solution x(t) as the Fourier series (5.94). 

In practice, one usually needs to include surprisingly few terms of the solution (5.94) 
to get a satisfactory approximation, as the following example illustrates. 20 

example 5.5 An Oscillator Driven by a Rectangular Pulse 

Consider a weakly damped oscillator that is being driven by the periodic rectan¬ 
gular pulses of Example 5.4 (Figure 5.22). Let the natural period of the oscillator 
be r 0 = 1, so that the natural frequency is co 0 — 2tc, and let the damping constant 


19 The n = 0, constant term needs separate consideration. It is easy to see that for a constant 
force / 0 the solution is x 0 — fo/cof This is actually exactly what you get if you just set n = 0 in 
(5.92) and (5.93). 

20 The solution contained in Equations (5.92) to (5.94) can be written more compactly if you 
don’t mind using complex notation. See Problem 5.51. 
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be ft = 0.2. Let the pulse last for a time At = 0.25 and have a height / max = 1. 
Calculate the first six Fourier coefficients A n for the long-term motion x (t) of the 
oscillator, assuming first that the drive period is the same as the natural period, 
r = r G = 1. Plot the resulting motion for several complete oscillations. Repeat 
these exercises for r = 1.5r 0 , 2.0 t o , and 2.5r 0 . 

Before we look at any of these exercises, it is worth thinking of a real system 
this problem might represent. One simple possibility is a mass hanging from the 
end of a spring, to which a professor is applying regularly spaced upward taps 
at intervals r. An even more familiar example would be a child in a swing, to 
whom a parent is giving regularly spaced pushes — though in this case we need 
to be careful to keep the amplitude small to justify the use of Hooke’s law. We 
are told to start by taking r = r G = 1; that is, the parent is pushing the child at 
exactly the natural frequency. 

The Fourier coefficients f n of the driving force were already calculated in 
Equations (5.86) and (5.87) of Example 5.4 (where they were called a n ). If we 
substitute these into (5.92) for the coefficients A n and put in the given numbers 
(including r = r 0 = 1), we find for the first six Fourier coefficients A 0 , • • •, A 5 : 

| 

Ai A 2 A 3 A 4 A 5 

63 1791 27 5 0 1 

(Since the numbers are rather small, I have quoted the values multiplied by 10 4 ; 
that is, A 0 = 63 x 10~ 4 , A l = 1791 x 10~ 4 , and so on.) Two things stand out 
about these numbers: First, after A h they get rapidly smaller, and for almost all 
purposes it would be an excellent approximation to ignore all but the first three 
terms in the Fourier series for x{t). Second, the coefficient A { is vastly bigger 

i 

1 than all the rest. This is easy to understand if you look at (5.92) for the coefficient 
Ap Since oo = co 0 (remember the parent is pushing the child at the swing’s 

I 

natural frequency) and n = 1, the first term in the denominator is exactly zero, 
the denominator is anomalously small, and A x is anomalously big compared 
to all the other coefficients. In other words, when the driving frequency is the 
same as the natural frequency, the n = 1 term in the Fourier series for x(t) is at 
resonance , and the oscillator responds especially strongly with frequency co 0 . 

Before we can plot x{t) as given by (5.94), we need to calculate the phase 
shifts 8 n using (5.93). This is easily done, though I shan’t waste space displaying 
the results. We cannot actually plot the infinite series (5.94); instead, we must 
pick some finite number of terms with which to approximate x(t). In the present 
case, it seems clear that three terms would be plenty, but to be on the safe side 
I’ll use six. Figure 5.24 shows x(t) as approximated by the sum of the first 
six terms in (5.94). At the scale shown, this approximate graph is completely 
indistinguishable from the exact result, which in turn is indistinguishable from 
a pure cosine 21 with frequency equal to the natural frequency of the oscillator. 
The strong response at the natural frequency is just what we would expect. For 


21 Actually it’s a pure sine, but this is really cos (cot — 8 { ) with 8^ = n/2 as we should have 
expected because we’re exactly on resonance. 










Section 5.8 Fourier Series Solution for the Driven Oscillator* 


201 



Figure 5.24 The motion of a linear oscillator, driven by periodic 
rectangular pulses, with the drive period r equal to the natural 
period r 0 of the oscillator (and hence co = co 0 ). The horizontal 
axis shows time in units of the natural period r 0 . As expected the 
motion is almost perfectly sinusoidal, with period equal to the 
natural period. 

instance, anyone who has pushed a child on a swing knows that the most efficient 
way to get the child swinging high is to administer regularly spaced pushes at 
intervals of the natural period — that is, r = r 0 — and that the swing will then 
oscillate vigorously at its natural frequency. 

A driving force with any other period r can be treated in exactly the same 
way. The Fourier coefficients A 0 , • • •, A 5 for all of the values of r requested 
above are shown in Table 5.1. 


Table 5.1 The first six Fourier coefficients A n for the motion x(t) 
of a linear oscillator driven by periodic rectangular pulses, for four 
different drive periods r = r Q , 1.5r 0 , 2.0r o , and 2.5r 0 . All values have 
been multiplied by 10 4 . 




^0 

A i 

A 2 

A 3 

a 4 

^5 

X 

K>° 

o 

ii 

63 

1791 

27 

5 

0 

1 

X 

ii 

Ln 

o H 

42 

145 

89 

18 

6 

2 

X 

= 2.0 r 0 

32 

82 

896 

40 

13 

6 

X 

= 2.5 r 0 

25 

59 

130 

97 

25 

11 


The entries in the four rows of this table deserve careful examination. The 
first row (t = r 0 ) shows the coefficients already discussed, the most prominent 
feature of which is that the n = 1 coefficient is far the largest, because it is 
exactly on resonance. In the next row (r = 1.5r 0 ), the n = 1 Fourier component 
has moved well away from resonance, and A x has dropped by a factor of 12 
or so. Some of the other coefficients have increased a bit, but the net effect is 
that the oscillator moves much less than when r = r 0 . This is clearly visible in 
Figures 5.25(a) and (b), which show x(t) (as approximated by the first six terms 
of its Fourier series) for these two values of the drive period. 



202 


Chapter 5 Oscillations 





Figure 5.25 The motion of a linear oscillator, driven by periodic rectangular pulses, showing four 
different values of the drive period r. (a) When the oscillator is driven at its natural frequency (r = r 0 ), 
the n = 1 term of the Fourier series is at resonance and the oscillator responds energetically, (b) When 
t = 1.5r 0 , the response is feeble, (c) When r = 2.0r o , the n = 2 term is at resonance and the response 
is strong again, (d) When r = 2.5r 0 , the response is again weak. 


The third row of Table 5.1 shows the Fourier coefficients for a drive period 
equal to twice the natural period — the parent is pushing just once every second 
swing of the child. Now, with r = 2.0t o , the drive frequency is half the natural 
frequency (co = \co Q ). This means that the n = 2 Fourier component, with fre¬ 
quency 2 co = co 0 , is exactly on resonance, and the coefficient A 2 is anomalously 
large. Once again we get a large response, as seen in Figure 5.25(c). 

Let us look a little closer at the case of Figure 5.25(c). It is, of course, a 
matter of experience that a perfectly satisfactory way to get a child swinging 
is to push once every two swings, although this naturally doesn’t get quite the 
result of pushing once for every swing. If you look carefully at Figure 5.25(c), 
you will notice that the swings alternate in size — the even swings are slightly 
bigger than the odd ones. This too is to be expected: Because the oscillator is 
damped, the second swing after each push is bound to be a little smaller than 
| the first. 

Finally, when r = 2.5r 0 , the n = 2 Fourier component is well past resonance 
and A 2 is much smaller again. On the other hand, the n = 3 component is 
approaching resonance so that A 3 is getting bigger. In fact, A 2 and A 3 are roughly 
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the same size, so that x(t) contains two dominant Fourier components and shows 
the somewhat kinky behavior seen in Figure 5.25(d). (Similar considerations 
apply to the case r = 1.5r 0 , where the coefficients A x and A 2 are both fairly 
large.) 


5.9 The RMS Displacement; Parseval’s Theorem* 

* The Parseval relation, which we introduce and apply in this section , is one of the most useful 
properties of Fourier series. Nevertheless, you could omit this section if pressed for time. 

In the last section we studied how the response of an oscillator varied with the 
frequency of the applied periodic driving force. We did this by solving for the motion 
x (t), using the method of Fourier series, for each of several interesting applied 
frequencies. It would be convenient if we could find a single number to measure the 
oscillator’s response and then just plot this number against the driving frequency (or 
driving period). In fact there are several ways to do this. Perhaps the most obvious 
thing to try would be the oscillator’s average displacement from equilibrium, (x ). (I am 
using angle brackets () to indicate a time average.) Unfortunately, since the oscillator 
spends as much time in any region of positive x as in the corresponding region of 
negative x, the average (x) is zero. 22 To get around this difficulty the most convenient 
quantity to use is the mean square displacement (x 2 ), and to give a quantity with the 
dimensions of length we usually discuss the root mean square or RMS displacement 

*rms = FF)- (5-95) 

The definition of the time average needs a little care. The usual practice is to define 
() as the average over one period r. Thus, 

i r ' 2 

(.x 2 ) = — / x 2 dt. (5.96) 

T J—r/2 

Because the motion is periodic, this is the same as the average over any integer number 
of periods, and hence also the average over any long time. (If this isn’t clear to you, 
see Problem 5.54.) 

To evaluate the average (x 2 ), we use the Fourier expansion (5.94) of x(t) 

00 

x(t) = A n cos (not — 8 n ). (5.97) 

n =0 

(In general, this series will contain sines as well as cosines, but in Example 5.5 the 
driving force contained only cosines, and for simplicity let us continue to assume this 


99 

Not quite true. The small constant term A 0 in (5.94) contributes a nonzero average (x) = A 0 , 
but this does not reflect any oscillation , which is what we are trying to characterize. 
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is the case. For the general case, see Problem 5.56.) Substituting for each of the factors 
x in (5.96), we get the appalling double sum 




A n cosincot 


8 n )A m cos(mcot — 8 m )dt. 


(5.98) 


Fortunately, this simplifies dramatically. It is fairly easy to show (Problem 5.55) that 
the integral is just 



cos (n cot — 8 n ) cos (mcot — 8 m ) dt = 


r if m = n = 0 

r/2 if m = n ^ 0 
0 if m n. 


(5.99) 


Thus, in the double sum (5.98), only those terms with m = n need to be retained, and 
we get the surprisingly simple result that 

oo 

<* 2 > = A 0 2 + i £>„ 2 . (5.100) 

n — 1 


This relation is called Parseval’s theorem . 23 It has many important theoretical uses, 
but for our purposes its main application is this: Since we know how to calculate 
the coefficients A n , Parseval’s theorem lets us find the response ( x 2 ) of our oscillator. 
Moreover, by dropping all but some modest finite number of terms in the sum (5.100), 
we get an excellent and easily calculated approximation for (x 2 ), as the following 
example illustrates. 


example 5.6 The RMS Displacement for a Driven Oscillator 

Consider again the oscillator of Example 5.5, driven by the periodic rectangular 
pulses of Example 5.4 (Figure 5.22). Find the RMS displacement = yj{x 2 ) 
as given by (5.100) for this oscillator. Using the same numerical values as 
before (r G = 1, = 0.2, / max — 1, Ar = 0.25) and approximating (5.100) by 

] its first six terms, make a plot of ^ rms as a function of the drive period r for 
| 0.25 < r < 5.5. 

I 

We have already done all the calculations needed to write down a formula 
I ^° r x rms — v ( x2 )- First, (x 2 ) is given by (5.100), where the Fourier coefficients 
A n are given by (5.92) as 

I 



_ fn _ 

~ n 2 a> 2 ) 2 + 4 P 2 n 2 co 2 


(5.101) 


23 Remember that we made the simplifying assumption that our Fourier series contained only 
cosine terms. In general, the sum in (5.100) must include contributions B 2 from the sine terms as 
well. See Problem 5.56. 
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Figure 5.26 The RMS displacement of a linear oscillator, driven by 
periodic rectangular pulses, as a function of the drive period r — 
calculated using the first six terms of the Parseval expression (5.100). 
The horizontal axis shows r in units of the natural period t 0 . When r 
is an integral multiple of r Q the response is especially strong. 


and the Fourier coefficients f n of the driving force are given by (5.87) and (5.86) 
as 


fn _ ^/max sin /tTwAt \ , | |or/; > ,| (5.102) ) 

it n \ r / 

while f 0 = / max Ar/r. Putting all these together gives the desired formula for 
x rms (which I’11 leave you to write down if you want to see it). 

If we now put in the given numbers, we are left with just one independent | 
variable, the period of the driving force r. (Remember that co = Trun¬ 

cating the infinite series (5.100) after six terms, we arrive at an expression that is 
easily evaluated with the appropriate software (or even a programmable calcula¬ 
tor) and plotted as shown in Figure 5.26. This graph shows clearly and succinctly 
what we found in the previous example. As we increase the drive period r, the | 
response of the oscillator varies dramatically. Each time r passes through an 
integer multiple of the natural period r Q (that is, r = nr Q ), the response exhibits j 
a sharp maximum, because the nth Fourier component is at resonance. On the j 
other hand, each successive peak is lower than its predecessor, since we elected 
to fix the width Ar and height / max of the pulses; thus as the drive period gets 
longer, the net effect of the force would be expected to get less. 


Principal Definitions and Equations of Chapter 5 

Hooke’s Law 


F = —kx <£=>■ U = | kx 2 


[Section 5.1] 
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Simple Harmonic Motion 

x = —arx <=> x(t) = A cos (cot — S), etc. [Section 5.2] 


Damped Oscillations 

If the oscillator is subject to a damping force —bv, then 

x + 2 fix + (OqX = 0 x(t) = Ae~^ cos {co x t — <$) [Eqs. (5.28) & (5.38)] 

where = 2 b/m, co Q = k/m , = ^Jco 2 _ ^ an( j solution given here is for 

“weak damping” < co 0 ). 


Driven Damped Oscillations and Resonance 

If the oscillator is also subject to a sinusoidal driving force F(t ) = mf 0 cos(cot), the 
long-term motion has the form 


x(t ) = Acos(cot — 5) 


[Eq. (5.66)] 


where 


A 


2 



(co 2 — co 2 ) 2 + 4 /3 2 co 2 


[Eq. (5.64)] 


and the phase shift 8 is given by (5.65). To this solution can be added a “transient” 
solution of the corresponding homogeneous equation, but this dies out as time passes. 
The long-term solution “resonates” (has a sharp maximum) when co is close to co 0 . 


Fourier Series 

If the driving force is not sinusoidal, but is still periodic, it can be built up as a Fourier 
series of sinusoidal terms, as in (5.90), and the resulting motion is the corresponding 
series of sinusoidal solutions, as in (5.94): 

(X) 

x(t ) = ^ A n cos (ncot — 8 n ). [Eq. (5.94)] 

n =0 


The RMS Displacement 

The root mean square displacement 



[Eqs. (5.95) & (5.96)] 
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is a good measure of the average response of the oscillator and is given by Parseval’s 
theorem as 


^rms 




00 


a? + !E a . 2 


n— 1 


[Eq. (5.100)] 


Problems for Chapter 5 

Stars indicate the approximate level of difficulty, from easiest (*) to most difficult (★★★). 


section 5.1 Hooke’s Law 

5.1 ★ A massless spring has unstretched length l 0 and force constant k. One end is now attached to the 
ceiling and a mass m is hung from the other. The equilibrium length of the spring is now f. (a) Write 
down the condition that determines l\. Suppose now the spring is stretched a further distance x beyond 
its new equilibrium length. Show that the net force (spring plus gravity) on the mass is F — —kx. That 
is, the net force obeys Hooke’s law, when x is the distance from the equilibrium position — a very useful 
result, which lets us treat a mass on a vertical spring just as if it were horizontal, (b) Prove the same 
result by showing that the net potential energy (spring plus gravity) has the form U (x) = const + \kx 2 . 

5.2 ★ The potential energy of two atoms in a molecule can sometimes be approximated by the Morse 
function, 


U(r) = A 


(,(«->■)/* _ i ) 2 _ i 


where r is the distance between the two atoms and A, R, and S are positive constants with S R. Sketch 
this function for 0 < r < oo. Find the equilibrium separation r 0 , at which U (r) is minimum. Now write 
r = r Q + x so that x is the displacement from equilibrium, and show that, for small displacements, U 
has the approximate form U = const + \kx 2 . That is, Hooke’s law applies. What is the force constant kl 

5.3 ★ Write down the potential energy U (0) of a simple pendulum (mass m, length /) in terms of the 
angle 0 between the pendulum and the vertical. (Choose the zero of U at the bottom.) Show that, for 
small angles, U has the Hooke’s law form U (0) = \k(j) 2 , in terms of the coordinate 0. What is kl 

5.4 ★★ An unusual pendulum is made by fixing a string to a horizontal cylinder of radius R , wrapping 
the string several times around the cylinder, and then tying a mass m to the loose end. In equilibrium 
the mass hangs a distance l 0 vertically below the edge of the cylinder. Find the potential energy if the 
pendulum has swung to an angle 0 from the vertical. Show that for small angles, it can be written in 
the Hooke’s law form U = \k<f> 2 . Comment on the value of k. 


section 5.2 Simple Harmonic Motion 

5.5 ★ In Section 5.2 we discussed four equivalent ways to represent simple harmonic motion in one 
dimension: 
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x(t) = C x e ia)t + C 2 e~ icot (I) 

= B x cos(cot) + B 2 sin(cot) (II) 
= A cos (cot - 8) (III) 

= Re Ce icot (IV) 


To make sure you understand all of these, show that they are equivalent by proving the following 
implications: I =>► II => III => IV => I. For each form, give an expression for the constants (C h C 2 , etc.) 
in terms of the constants of the previous form. 

5.6 ★ A mass on the end of a spring is oscillating with angular frequency co. At t = 0, its position is 
x 0 > 0 and I give it a kick so that it moves back toward the origin and executes simple harmonic motion 
with amplitude 2x 0 . Find its position as a function of time in the form (III) of Problem 5.5. 

5.7 ★ (a) Solve for the coefficients B x and B 2 of the form (II) of Problem 5.5 in terms of the initial 
position jc 0 and velocity v Q at t = 0. (b) If the oscillator’s mass is m = 0.5 kg and the force constant 
is k = 50 N/m, what is the angular frequency col If x 0 = 3.0 m and v 0 = 50 m/s, what are B x and B 2 1 
Sketch x(t) for a couple of cycles, (c) What are the earliest times at which x = 0 and at which x = 0? 

5.8 ★ (a) If a mass m = 0.2 kg is tied to one end of a spring whose force constant k = 80 N/m and 
whose other end is held fixed, what are the angular frequency co , the frequency /, and the period r of its 
oscillations? (b) If the initial position and velocity are x o = 0 and v Q = 40 m/s, what are the constants 
A and 8 in the expression x(t) = A cos (cot — 8)1 

5.9 ★ The maximum displacement of a mass oscillating about its equilibrium position is 0.2 m, and its 
maximum speed is 1.2 m/s. What is the period r of its oscillations? 

5.10 ★ The force on a mass m at position x on the x axis is F = —F Q sinhax, where F 0 and a are 
constants. Find the potential energy U (v), and give an approximation for U(x) suitable for small 
oscillations. What is the angular frequency of such oscillations? 

5.11 ★ You are told that, at the known positions x x and jc 2 , an oscillating mass m has speeds v x and v 2 . 
What are the amplitude and the angular frequency of the oscillations? 

5.12 ★★ Consider a simple harmonic oscillator with period r. Let (/) denote the average value of any 
variable f(t), averaged over one complete cycle: 

(/> = - f f(t)dt. (5.103) 

t Jo 

Prove that (T) — (U) = \E where E is the total energy of the oscillator. [Hint: Start by proving the 
more general, and extremely useful, results that (sin 2 (cot — 8)) = (cos 2 (cot — <$)) = £. Explain why 
these two results are almost obvious, then prove them by using trig identities to rewrite sin 2 6 and 
cos 2 0 in terms of cos (20).] 

5.13 ★★ The potential energy of a one-dimensional mass m at a distance r from the origin is 

for 0 < r < oo, with £/ 0 , R, and X all positive constants. Find the equilibrium position r G . Let x be the 
distance from equilibrium and show that, for small x, the PE has the form U = const + \kx 2 . What is 
the angular frequency of small oscillations? 
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Figure 5.27 Problem 5.18 


section 5.3 Two-Dimensional Oscillators 

5.14 ★ Consider a particle in two dimensions, subject to a restoring force of the form (5.21). (The two 
constants k x and k y may or may not be equal; if they are, the oscillator is isotropic.) Prove that its 
potential energy is 

U = \{k x x 2 + k y y 2 ). (5.104) 

5.15 ★ The general solution for a two-dimensional isotropic oscillator is given by (5.19). Show that 
by changing the origin of time you can cast this in the simpler form (5.20) with 8 = 8 y — 8 X . [Hint: A 
change of origin of time is a change of variables from t to t' = t + t 0 . Make this change and choose 
the constant t 0 appropriately, then rename t' to be t .] 

5.16 ★ Consider a two-dimensional isotropic oscillator moving according to Equation (5.20). Show 
that if the relative phase is 8 = jt/2, the particle moves in an ellipse with semimajor and semiminor 
axes A x and A y . 

5.17 ★★ Consider the two-dimensional anisotropic oscillator with motion given by Equation (5.23). 
(a) Prove that if the ratio of frequencies is rational (that is, co x /a) y = p/q where p and q are integers) 
then the motion is periodic. What is the period? (b) Prove that if the same ratio is irrational, the motion 
never repeats itself. 

5.18 ★★★ The mass shown from above in Figure 5.27 is resting on a frictionless horizontal table. Each 
of the two identical springs has force constant k and unstretched length / 0 . At equilibrium the mass rests 
at the origin, and the distances a are not necessarily equal to / Q . (That is, the springs may already be 
stretched or compressed.) Show that when the mass moves to a position (jc, y), with x and y small, the 
potential energy has the form (5.104) (Problem 5.14) for an anisotropic oscillator. Show that if a < l Q 
the equilibrium at the origin is unstable and explain why. 

5.19 ★★★ Consider the mass attached to four identical springs, as shown in Figure 5.7(b). Each spring 
has force constant k and unstretched length / 0 , and the length of each spring when the mass is at its 
equilibrium at the origin is a (not necessarily the same as / G ). When the mass is displaced a small 
distance to the point (x, y), show that its potential energy has the form \k’r 2 appropriate to an isotropic 
harmonic oscillator. What is the constant k r in terms of k? Give an expression for the corresponding 
force. 


section 5.4 Damped Oscillations 

5.20 ★ Verify that the decay parameter — y ft 2 — go 2 for an overdamped oscillator (/3 > co 0 ) de¬ 
creases with increasing /3. Sketch its behavior for co Q < ft < oo. 

5.21 ★ Verify that the function (5.43), x(t) = te~^, is indeed a second solution of the equation of 
motion (5.28) for a critically damped oscillator (ft = co Q ). 

5.22 ★ (a) Consider a cart on a spring which is critically damped. At time t — 0, it is sitting at its 
equilibrium position and is kicked in the positive direction with velocity v Q . Find its position x(t ) for 
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all subsequent times and sketch your answer, (b) Do the same for the case that it is released from rest at 
position x = x 0 . In this latter case, how far is the cart from equilibrium after a time equal to x 0 = 2jt /co Q , 
the period in the absence of any damping? 

5.23 ★ A damped oscillator satisfies the equation (5.24), where F dmp = — bx is the damping force. Find 
the rate of change of the energy E — \mx 2 + \kx 2 (by straightforward differentiation), and, with the 
help of (5.24), show that dE/dt is (minus) the rate at which energy is dissipated by F dmp . 

5.24 ★ In our discussion of critical damping (P = &> 0 ), the second solution (5.43) was rather pulled 
out of a hat. One can arrive at it in a reasonably systematic way by looking at the solutions for 
P < co 0 and carefully letting ft -> co 0 , as follows: For ft < co 0 , we can write the two solutions as 
x x (t) — e~^ cos(o) 1 t) and x 2 (t) = e~^ sin (co x t ).'Show that as ft -> o> 0 , the first of these approaches 
the first solution for critical damping, x x (t) — e~^. Unfortunately, as ft co 0 , the second of them 
goes to zero. (Check this.) However, as long as ft ^ co 09 you can divide x 2 (t) by oo x and you will still 
have a perfectly good second solution. Show that as ft —> co Q , this new second solution approaches the 
advertised te~^. 

5.25 ★★ Consider a damped oscillator with p < co 0 . There is a little difficulty defining the “period” x x 
since the motion (5.38) is not periodic. However, a definition that makes sense is that x x is the time 
between successive maxima of jt(t). (a) Make a sketch of x(t) against t and indicate this definition of 
r on your graph. Show that x x — 2 7t/co x . (b) Show that an equivalent definition is that x x is twice the 
time between successive zeros of x(t). Show this one on your sketch, (c) If P = co 0 / 2, by what factor 
does the amplitude shrink in one period? 

5.26 ★★ An undamped oscillator has period x 0 = 1.000 s, but I now add a little damping so that its 
period changes to x x = 1.001 s. What is the damping factor pi By what factor will the amplitude of 
oscillation decrease after 10 cycles? Which effect of damping would be more noticeable, the change 
of period or the decrease of the amplitude? 

5.27 ★★ As the damping on an oscillator is increased there comes a point when the name “oscillator” 
seems barely appropriate, (a) To illustrate this, prove that a critically damped oscillator can never pass 
through the origin x = 0 more than once, (b) Prove the same for an overdamped oscillator. 

5.28 ★★ A massless spring is hanging vertically and unloaded, from the ceiling. A mass is attached to 
the bottom end and released. How close to its final resting position is the mass after 1 second, given that 
it finally comes to rest 0.5 meters below the point of release and that the motion is critically damped? 

5.29 ★★ An undamped oscillator has period x 0 = 1 second. When weak damping is added, it is found 
that the amplitude of oscillation drops by 50% in one period x x . (The period of the damped oscillations 
is defined as the time between successive maxima, x x = 2n/co x . See Problem 5.25.) How big is P 
compared to oo 0 l What is x{l 

5.30 ★★ The position x(t) of an overdamped oscillator is given by (5.40). (a) Find the constants C x 
and C 2 in terms of the initial position x 0 and velocity v 0 . (b) Sketch the behavior of x(t) for the two 
cases that v Q = 0 and that x o = 0. (c) To illustrate again how mathematics is sometimes cleverer than 
we (and check your answer), show that if you let P -> 0, your solution for x(t) in part (a) approaches 
the correct solution for undamped motion. 

5.31 ★★ [Computer] Consider a cart on a spring with natural frequency co 0 = 2n, which is released 
from rest at x 0 = 1 and t = 0. Using appropriate graphing software, plot the position x(t) for 0 < t < 2 
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and for damping constants f3 =0, 1, 2, 4, 6, 27t, 10, and 20. [Remember that x(t) is given by different 
formulas for < co 0 , f$ = co 0 , and > co 0 .\ 

5.32 ★★ [Computer] Consider an underdamped oscillator (such as a mass on the end of a spring) that 
is released from rest at position x G at time / = 0. (a) Find the position x(t) at later times in the form 

x(t) = cos(co x t) + B 2 sin(o)\t)]. 

That is, find B { and B 2 in terms of x 0 . (b) Now show that if you let approach the critical value co 0 , 
your solution automatically yields the critical solution, (c) Using appropriate graphing software, plot 
the solution for 0 < t < 20, with x 0 = 1, co 0 = 1, and = 0, 0.02, 0.1, 0.3, and 1. 

section 5.5 Driven Damped Oscillations 

5.33 ★ The solution for x(t) for a driven, underdamped oscillator is most conveniently found in the 
form (5.69). Solve that equation and the corresponding expression for i, to give the coefficients B x 
and B 2 in terms of A, <5, and the initial position and velocity x 0 and v 0 . Verify the expressions given in 
(5.70). 

5.34 ★ Suppose that you have found a particular solution x p (t) of the inhomogeneous equation (5.48) 
for a driven damped oscillator, so that Dx p = / in the operator notation of (5.49). Suppose also 
that x(t ) is any other solution, so that Dx = /. Prove that the difference x — x p must satisfy the 
corresponding homogeneous equation, D (x — x p ) = 0. This is an alternative proof that any solution x of 
the inhomogeneous equation can be written as the sum of your particular solution plus a homogeneous 
solution; that is, x = x p + x h . 

5.35 ★★ This problem is to refresh your memory about some properties of complex numbers needed 
at several points in this chapter, but especially in deriving the resonance formula (5.64). (a) Prove that 
any complex number z = x + iy (with x and y real) can be written as z = re 16 where r and 6 are the 
polar coordinates of z in the complex plane. (Remember Euler’s formula.) (b) Prove that the absolute 
value of z, defined as |z| = r, is also given by |z| 2 = zz*, where z* denotes the complex conjugate of z, 
defined as z* = x — iy. (c) Prove that z* = re ~ ld . (d) Prove that (zw)* = z*w* and that (1/z)* = 1/z*. 
(e) Deduce that if z = a/(b + zc), with a , b , and c real, then |z| 2 = a 2 /(b 2 + c 2 ). 

5.36 ★★ [Computer] Repeat the calculations of Example 5.3 (page 185) with all the same parameters, 
but with the initial conditions x 0 = 2 and v 0 = 0. Plot x(t) for 0 < t < 4 and compare with the plot of 
Example 5.3. Explain the similarities and differences. 

5.37 ★★ [Computer] Repeat the calculations of Example 5.3 (page 185) but with the following param¬ 
eters 


co = 2 tt, co 0 = 0.25 co, = 0.2 co 0 , f 0 = 1000 

and with the initial conditions x 0 = 0 and v Q = 0. Plot x(t) for 0 < t < 12 and compare with the plot 
of Example 5.3. Explain the similarities and differences. (It will help your explanation if you plot the 
homogeneous solution as well as the complete solution — homogeneous plus particular.) 

5.38 ★★ [Computer] Repeat the calculations of Example 5.3 (page 185) but take the parameters of the 
system to be co = co 0 = 1, /3 = 0.1, and f Q = 0.4, with the initial conditions x 0 = 0 and v Q = 6 (all in 
some apppropriate units). Find A and 8 , and then B x and B 2 , and make a plot of x(t) for the first ten or 
so periods. 
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5.39 ★★ [Computer] To get some practice at solving differential equations numerically, repeat the 
calculations of Example 5.3 (page 185), but instead of finding all the various coefficients just use 
appropriate software (for example, the NDSolve command of Mathematica) to solve the differential 
equation (5.48) with the boundary conditions x 0 = v 0 = 0. Make sure your graph agrees with Figure 
5.15. 

section 5.6 Resonance 

5.40 ★ Consider a damped oscillator, with fixed natural frequency co Q and fixed damping constant f$ 
(not too large), that is driven by a sinusoidal force with variable frequency co. Show that the amplitude 

of the response, as given by (5.71) is maximum when co = jco 2 — 2/3 2 . (Note that so long as the 
resonance is narrow this implies co & co Q .) 

5.41 ★ We know that if the driving frequency co is varied, the maximum response (A 2 ) of a driven 
damped oscillator occurs at co ^oo 0 (if the natural frequency is co 0 and the damping constant /3 
co Q ). Show that A 2 is equal to half its maximum value when co & co Q =t /?, so that the full width at 
half maximum is just 2/3. [Hint: Be careful with your approximations. For instance, it’s fine to say 
co + co 0 as 2 co 0 , but you certainly mustn’t say co — co 0 ~ ().] 

5.42 ★ A large Foucault pendulum such as hangs in many science museums can swing for many hours 
before it damps out. Taking the decay time to be about 8 hours and the length to be 30 meters, find the 
quality factor Q. 

5.43 ★★ When a car drives along a “washboard” road, the regular bumps cause the wheels to oscillate 
on the springs. (What actually oscillates is each axle assembly, comprising the axle and its two wheels.) 
Find the speed of my car at which this oscillation resonates, given the following information: (a) When 
four 80-kg men climb into my car, the body sinks by a couple of centimeters. Use this to estimate the 
spring constant k of each of the four springs, (b) If an axle assembly (axle plus two wheels) has total 
mass 50 kg, what is the natural frequency of the assembly oscillating on its two springs? (c) If the 
bumps on a road are 80 cm apart, at about what speed would these oscillations go into resonance? 

5.44 ★★ Another interpretation of the Q of a resonance comes from the following: Consider the motion 
of a driven damped oscillator after any transients have died out, and suppose that it is being driven close 
to resonance, so you can set co = co 0 . (a) Show that the oscillator’s total energy (kinetic plus potential) 
is E = \mco 2 A 2 . (b) Show that the energy AF dis dissipated during one cycle by the damping force 
F dmp is 2jtmPcoA 2 . (Remember that the rate at which a force does work is Fv.) (c) Hence show that 
Q is 2jz times the ratio E / A£ dis . 

5.45 ★★★ Consider a damped oscillator, with natural frequency co Q and damping constant both fixed, 
that is driven by a force F(t ) = F 0 cos (cot), (a) Find the rate P(t) at which F(t) does work and show 
that the average rate (P) over any number of complete cycles is mpco 2 A 2 . (b) Verify that this is the 
same as the average rate at which energy is lost to the resistive force, (c) Show that as co is varied {P) 
is maximum when co = co Q \ that is, the resonance of the power occurs at co = co 0 (exactly). 

section 5.7 Fourier Series * 

5.46 ★ The constant term a 0 in a Fourier series is a bit of a nuisance, always requiring slightly special 
treatment. At least it has a rather simple interpretation: Show that if f(t) has the standard Fourier series 
(5.82), then a 0 is equal to the average (/) of f(t) taken over one complete cycle. 
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5.47 ★★ In order to prove the crucial formulas (5.83)11(5.85) for the Fourier coefficients a n and b n , you 
must first prove the following: 



cos {ncot) cos (moot) dt — 


r/2 if m = n ^ 0 
0 if m n. 


(5.105) 


(This integral is obviously r if m = n = 0.) There is an identical result with all cosines replaced by 
sines, and finally 

r */2 

I cos (ncot) sin {moot) dt = 0 for all integers n and m, (5.106) 

J- r/2 

where as usual co — Iti/x. Prove these. [Hint: Use trig identities to replace cos(0) cos(</>) by terms like 
co$(0 + 0) and so on.] 

5.48 ★★ Use the results (5.105) and (5.106) to prove the formulas (5.83)—(5.85) for the Fourier coeffi¬ 
cients a n and b n . [Hint: Multiply both sides of the Fourier expansion (5.82) by cos {moot) or sin (moot) 
and then integrate from —r/2 to r/2.] 

5.49 ★★★ [Computer] Find the Fourier coefficients a n and b n for the function shown in Figure 5.28(a). 
Make a plot similar to Figure 5.23, comparing the function itself with the first couple of terms in the 
Fourier series, and another for the first six or so terms. Take / max = 1. 




Figure 5.28 (a) Problem 5.49. (b) Problem 5.50 


5.50 ★★★ [Computer] Find the Fourier coefficients a n and b n for the function shown in Figure 5.28(b). 
Make a plot similar to Figure 5.23, comparing the function itself with the sum of the first couple of 
terms in the Fourier series, and another for the first 10 or so terms. Take / max = 1. 

section 5.8 Fourier Series Solution for the Driven Oscillator* 

5.51 ★★ You can make the Fourier series solution for a periodically driven oscillator a bit tidier if you 
don’t mind using complex numbers. Obviously the periodic force of Equation (5.90) can be written as 
/ = Re(g), where the complex function g is 

oo 

g(o = £/y' na ". 

n= 0 

Show that the real solution for the oscillator’s motion can likewise be written as x = Re(z), where 

oo 

z(t) = J2 C n e inM 

n =0 
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and 

^ _ fYl _ 

n co 2 — n 2 co 2 + lifinod 

This solution avoids our having to worry about the real amplitude A n and phase shift 8 n separately. (Of 
course A n and 8 n are hidden inside the complex number C n .) 

5.52 ★★★ [Computer] Repeat all the calculations and plots of Example 5.5 (page 199) with all the same 
parameters except that =0.1. Compare your results with those of the example. 

5.53 ★★★ [Computer] An oscillator is driven by the periodic force of Problem 5.49 [Figure 5.28(a)], 
which has period r = 2. (a) Find the long-term motion x (f), assuming the following parameters: natural 
period r 0 = 2 (that is, co 0 = 7r), damping parameter /3 =0.1, and maximum drive strength / max = 1. 
Find the coefficients in the Fourier series for x(t) and plot the sum of the first four terms in the series 
for 0 < t < 6. (b) Repeat, except with natural period equal to 3. 

section 5.9 The RMS Displacement; Parseval’s Theorem* 

5.54 ★ Let f(t) be a periodic function with period r. Explain clearly why the average of / over one 
period is not necessarily the same as the average over some other time interval. Explain why, on the 
other hand, the average over a long time T approaches the average over one period, as T -» oo. 

5.55 ★★ To prove the Parseval relation (5.100), one must first prove the result (5.99) for the integral of 
a product of cosines. Prove this result, and then use it to prove the Parseval relation. 

5.56 ★★ The Parseval relation as stated in (5.100) applies to a function whose Fourier series happens 
to contain only cosines. Write down the relation and prove it for a function 

oo 

x(0 = [ A n cos (ncot — 8 n ) -F B n sin(n6L>? — <$„)]. 

n= 0 

5.57 ★★ [Computer] Repeat the calculations that led to Figure 5.26, using all the same parameters 
except taking /3 =0.1. Plot your results and compare your plot with Figure 5.26. 
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Calculus of Variations 


In many problems one needs to use non-Cartesian coordinates. Roughly speaking there 
are two classes of such problems. First, certain symmetries make it most advantageous 
to use special coordinates: Problems with spherical symmetry call out for the use of 
spherical polar coordinates; similarly, problems with axial symmetry are best treated 
in cylindrical polar coordinates. Second, when particles are constrained in some way, 
it is usually best to choose an appropriate, and usually non-Cartesian, coordinate 
system. For example, an object that is constrained to move on the surface of a sphere 
is probably best treated using spherical polar coordinates; if a bead slides on a curved 
wire, the best choice of coordinate may be just the distance along the curving wire 
from some convenient origin. 

Unfortunately, as we have seen, the expressions for the components of the accel¬ 
eration in non-Cartesian coordinates are quite messy, and the situation gets rapidly 
worse as we move on to more complicated systems. This makes Newton’s second 
law difficult to use in non-Cartesian coordinates. We need an alternative (though ul¬ 
timately equivalent) equation of motion that works equally well in any coordinates, 
and the required alternative is provided by Lagrange’s equations. 

The best way to prove — and to understand the great flexibility of — Lagrange’s 
equations is to use a “variational principle.” Variational principles are important in 
many areas of mathematics and physics. It has proved possible to formulate almost 
every branch of physics — classical mechanics, quantum mechanics, optics, electro¬ 
magnetism, and so on — in variational terms. To the beginning student, accustomed 
to Newton’s laws, a reformulation of classical mechanics in terms of a variational 
principle does not necessarily seem like an improvement. But because they allow a 
similar formulation of so many different subjects, variational methods have given a 
unity to physics and have played a crucial role in the recent history of physical theory. 
For this reason, I would like to introduce variational methods in a reasonably general 
setting. Therefore this short chapter is a brief introduction to variational problems in 
general. In the next chapter I shall apply what we learn here to establish the Lagrangian 
formulation of mechanics. If you are already familiar with the “calculus of variations” 
you could skip straight to Chapter 7. 
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6.1 Two Examples 


The calculus of variations involves finding the minimum or maximum of a quantity 
that is expressible as an integral. To see how this can arise, I would like to start with 
two simple, concrete examples. 


The Shortest Path between Two Points 


My first example is this problem: Given two points in a plane, what is the shortest 
path between them? While you certainly know the answer — a straight line — you 
probably have not seen a proof, unless you have studied the calculus of variations. 
The problem is illustrated in Figure 6.1, which shows the two given points, (x h jj) 
and ( x 2 , y 2 ), and a path, y = y (x), joining them. Our task is to find the path y(x) that 
has the shortest length and to show that it is in fact a straight line. 

The length of a short segment of the path is ds = yfdx 2 + dy 2 , which, since 


dy 


—dx = y'(x) dx , 
dx 


we can rewrite as 


ds = yjdx 2 + dy 2 = yjl + y'(x) 2 dx. (6.1) 

Thus the total length of the path between points 1 and 2 is 

/ 2 rx 2 

ds — J y/l + y'{x) 2 dx. (6.2) 

This equation puts our problem in mathematical form: The unknown is the function 
y = y(x) that defines the path between points 1 and 2. The problem is to find the 
function y (x) for which the integral (6.2) is a minimum. It is interesting to contrast this 
with the standard minimization problem of elementary calculus, where the unknown 



Figure 6.1 A path joining the two points 1 and 2. The length of 

the short segment is ds = y/dx 2 + dy 2 , and the total length of 

2 

the path is L = f x ds. 
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is the value of a variable x at which a known function /(x) is a minimum. Obviously 
our new problem is one stage more complicated than this old one. 

Before we set up the machinery to solve this new problem, let’s consider another 
example. 


Fermat’s Principle 


A similar problem is to find the path that light will follow between two points. If the 
refractive index of the medium is constant, then the path is, of course, a straight line, 
but if the refractive index varies, or if we interpose a mirror or lens, the path is not so 
obvious. The French mathematician Fermat (1601-1665) discovered that the required 
path is the path for which the time of travel of the light is minimum. We can illustrate 
Fermat’s principle using Figure 6.1. The time for light to travel a short distance ds 
is ds/v where v denotes the speed of light in the medium, v = c/n where n is the 
refractive index. Thus Fermat’s principle says that the correct path between points 1 
and 2 is the path for which the time 


(time of travel) — 



1 

c 



n ds 


is a minimum. If n is constant, then it can be taken outside the integral and the problem 
reduces to finding the shortest path between points 1 and 2 (and the answer is, of 
course, a straight line). In general, the refractive index can vary, n = rc(x, y), and our 
problem is to find the path y(x) for which the integral 



n(x, y) ds 



n(x , y)y/\ + y'(x) 2 dx 


(6.3) 


is minimum. [In writing the last expression, I substituted (6.1) for ds.] 

The integral that has to be minimized in connection with Fermat’s principle is 
very similar to the integral (6.2) giving the length of a path; it is just a little more 
complicated, since the factor n(x, y) introduces an extra dependence on x and y. 
Similar integrals arise in many other problems. Sometimes we want the path for 
which an integral is a maximum , and sometimes we are interested in both maxima 
and minima. To get some idea of the possibilities, it is helpful to think again about the 
problem of finding maxima and minima of functions in elementary calculus. There we 
know that the necesary condition for a maximum or minimum of a function f(x) is 
that its derivative vanish, df/dx — 0. Unfortunately, this condition is not quite enough 
to guarantee a maximum or minimum. As you certainly recall from introductory 
calculus, there are essentially three possibilities, as illustrated in Figure 6.2. A point 
x 0 where df/dx is zero may be a maximum or a minimum or, if d 2 f/dx 2 is also zero, 
it may be neither , as indicated in Figure 6.2(c). When df/dx = 0 at a point x 0 , but 
we don’t know which of the three possibilities obtains, we say that x 0 is a stationary 
point of the function /(x), since an infinitesimal displacement of x from x G leaves 
/(x) unchanged (because the slope is zero). 

The situation for the problems of this chapter is very similar. The method I shall 
describe in the next section actually finds the path that makes an integral like (6.2) or 
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(a) (b) 



Figure 6.2 If df/dx = 0 at x 0 , there are three possibilities: (a) If the second 
derivative is positive, then f(x ) has a minimum at x 0 . (b) If the second 
derivative is negative, then f(x) has a maximum, (c) If the second derivative 
is zero, then there may be a minimum, a maximum, or neither (as shown). 


(6.3) stationary, in the sense that an infinitesimal variation of the path from its correct 
course doesn’t change the value of the integral concerned. If you need to know that the 
integral is definitely minimum (or definitely maximum, or perhaps neither), you have 
to check this separately. Incidentally, we are now ready to explain the name of this 
chapter: Since our concern is how infinitesimal variations of a path change an integral, 
the subject is called the calculus of variations. For the same reason, the methods we 
shall develop are called variational methods, and a principle like Fermat’s principle 
is a variational principle. 


6.2 The Euler-Lagrange Equation 


The two examples of the last section illustrate the general form of the so-called 
variational problem. We have an integral of the form 



where y(x) is an as-yet unknown curve joining two points (x h yj) and (x 2 , y?) as in 
Figure 6.1; that is, 


y(x i) = y x and y(x 2 ) = y 2 . (6.5) 

Among all the possible curves satisfying (6.5) (that is, joining the points 1 and 2), 
we have to find the one that makes the integral S a minimum (or maximum or at least 
stationary). To be definite, I shall suppose that we wish to find a minimum. Notice that 
the function / in (6.4) is a function of three variables / = f(y, y', x), but because 
the integral follows the path y = y(x) the integrand f[y(x), y'(x), x] is actually a 
function of just the one variable x. 
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Figure 6.3 The path y = y(x) between points 1 and 2 is the “right” 
path, the one for which the integral S of (6.4) is a minimum. Any other 
path Y (jc) is “wrong,” in that it gives a larger value for S. 


Let us denote the correct solution to our problem by y = y(x). Then the integral 
S in (6.4) evaluated for y = y (jc) is less than for any neighboring curve y = Y(x), as 
sketched in Figure 6.3. It is convenient to write the “wrong” curve Y(x) as 

Y(x) = j(x) + rj(x) (6.6) 

where rj(x) (Greek “eta”) is just the difference between the wrong Y (jc) and the right 
y(jt). Since Y(x) must pass through the endpoints 1 and 2, rj(x) must satisfy 

V(x l ) = r)(x 2 )=0. (6.7) 

There are infinitely many choices for the difference rj(x)\ for example, we could 
choose rj = (x — x x )(x 2 — x) or rj(x ) = sin[7r(x — x x )/{x 2 — x x )]. 

The integral S taken along the wrong curve Y(x) must be larger than that along 
the right curve y(v), no matter how close the former is to the latter. To express this 
requirement, I shall introduce a parameter a and redefine T(x) to be 

Y(x) = y(x) + arj(x). (6.8) 

The integral S taken along the curve F(x) now depends on the parameter a, so I shall 
call it S(a). The right curve y(x) is obtained from (6.8) by setting a = 0. Thus the 
requirement that S is minimum for the right curve y(x) implies that S'(a) is a minimum 
at a = 0. With this result, we have converted our problem to the traditional problem 
from elementary calculus of making sure that an ordinary function [namely S(a )] has 
a minimum at a specified point (a =0). To ensure this, we must just check that the 
derivative dS/da is zero when a = 0. 

If we write out the integral S(a) in detail, it looks like this: 

r* 2 

S(o0= / f(Y,Y\x)dx 

J X\ 

px 2 

= / f(y + ccri, / + ar}',x)dx. 

J Xi 


(6.9) 
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To differentiate (6.9) with respect to a, we note that a appears in the integrand /, so 
we need to evaluate df/da. Since a appears in two of the arguments of /, this gives 
two terms, namely (using the chain rule) 

df {y + art, y' + ari', x) df ,df 

-= r) -b r] —, 

da dy By' 


and for dS/da (which has to be zero) 
dS f x 2 df 


da 


j 

J X\ 


/.'(■ 


dx — I I rj— + r I dx = 0. 
da Jr , \ dy dy' 


( 6 . 10 ) 


This condition must be true for any rj(x) satisfying (6.7); that is, for any choice of the 
“wrong” path F(x) = y(x) + arj(x). 

To take advantage of the condition (6.10), we need to rewrite the second term on 
the right using integration by parts 1 (remember that r}' means drj/dx): 


px.2 

/ V 

J X\ 


(x)-f-dx = 
dy' 


rj{x) 


dy' J 


x 2 


r X2 d 

- / T,(X) 

Xj J X\ 


df\ 


dx 


dx \dy'J 


Because of the condition (6.7), the first term on the right (the “endpoint term”) is zero. 
Thus 2 


px 2 

/ V 

J X\ 


(x)^-dx = 
dy' 


f. 


X2 , , d (df - , 

rj(x) — — I dx. 

x, dx \dy' 


( 6 . 11 ) 


Substituting this identity into (6.10), we find that 

■x 2 


r m(f - 

Jx x \dy dx dy'J 


( 6 . 12 ) 


This condition must be satisfied for any choice of the function r\{x). Therefore, as I 
shall argue in a moment, the factor in large parentheses must be zero: 


0 (Buler-Lagraoge: 


for all x (in the relevant interval < x < x 2 ). This is the so-called Euler-Lagrange 
equation (named for the Swiss mathematician Leonhard Euler, 1707-1783, and the 
Italian-French physicist and mathematician Joseph Lagrange, 1736-1813), which lets 


1 If you are used to thinking of integration by parts in the form / v du = [uv] — / udv, then you 
will find it helpful to recognize that another way to say the same thing is: / u'vdx = [uv] — f uv' dx. 
In words: In the integral j u'v dx, you can move the prime from the u to the v if you change the 
sign and add the endpoint contribution [uv\ 

2 This is the simple form in which integration by parts often appears in physics: Provided the 
endpoint term [uv] is zero (as often happens), integration by parts lets you move the differentiation 
from the u to the v as long as you change the sign; that is, / u'v dx = —f uv'dx. 
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us find the path for which the integral S is stationary. Before I illustrate its use, I need 
to discuss the step from (6.12) to (6.13), which is by no means obvious. 

Equation (6.12) has the form f rj(x)g(x) dx — 0.1 would certainly not claim that 
this condition alone implies that g(;t) = 0 for all x. However, (6.12) holds for any 
choice of the function r](x ), and if J rj(x)g(x)dx = 0 for any rj(x ), then we can 
conclude that g(x) =0 for all x. To prove this, we must assume that all functions 
concerned are continuous, but, as physicists, we would take for granted that this is 
the case. 3 Now, to prove the assertion, let us assume the contrary, that g(x) is nonzero 
in some interval between x { and x 2 . Then choose a function rj(x) that has the same 
sign as g(x) (that is, rj is positive where g is positive and rj is negative where g is 
negative). Then the integrand is continuous, satisfies rj(x)g(x) > 0, and is nonzero 
at least in some interval. Under these conditions / rj(x)g(x) dx cannot be zero. This 
contradiction implies that g(x) is zero for all x. 

This completes the proof of the Euler-Lagrange equation. The procedure for using 
it is this: (1) Set up the problem so that the quantity whose stationary path you seek 
is expressed as an integral in the standard form 

s= f f[y(x)>y'(x),x]dx, (6.14) 

J X\ 

where f[y(x), y'(x ), x] is the function appropriate to your problem. (2) Write down 
the Euler-Lagrange equation (6.13) in terms of the function f[y(x ), y'(x),x]. (3) 
Finally, solve (if possible) the differential equation (6.13) for the function y(x) that 
defines the required stationary path. I shall illustrate this procedure with a couple of 
examples in the next section. 


6.3 Applications of the Euler-Lagrange Equation 


Let us start with the problem that began this chapter, finding the shortest path between 
two points in a plane. 

example 6.1 Shortest Path between Two Points 

We saw that the length of a path between points 1 and 2 is given by the integral 
(6.2) as 

/ 2 px 2 _ j 

ds = J yj 1 + y' 2 dx. 

This has the standard form (6.14), with the function / given by 

f(y,y',x) = (l + y' 2 ) l/2 . (6.15) i 


O 

The claimed result is clearly false if discontinuous functions are admitted. For instance, if we 
made g(x) nonzero at just one point, then J r](x)g(x) dx would still be zero. 
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To use the Euler-Lagrange equation (6.13), we must evaluate the two partial 
derivatives concerned: 


a/ n , 3/ y 

ay dy' (1 + y 2 )V 2 


Since df/dy = 0, (6.13) implies simply that 

d df 


dx dy' 


= 0 . 


(6.16) 


| In other words, df/dy' is a constant, C. According to (6.16), this implies that 
! y a — C 2 (l + y' 2 ), 

I 

or, with a little rearrangement, y a — constant. This implies that y\x) is a 
constant, which we could call m. Integrating the equation y'(x) = m, we find 
that y(x) = mx + b , and we have proved that the shortest path between two 
points is a straight line! 


A Note on Variables 

So far we have considered problems with two variables, which we have called x and 
y. Of these, x has been the independent variable, and y the dependent, through the 
relation y = y(x). Unfortunately, we are frequently forced — by convenience or tra¬ 
dition — to name the variables differently. For example, in a simple one-dimensional 
mechanics problem, the independent variable is the time t and the dependent variable 
is the position x = x(t). This means you will have to get used to seeing the Euler- 
Lagrange equation with the variables x and y replaced by an assortment of other 
variables, such as t and x. In the next example, the two variables are x and y, but the 
independent variable is y, and the roles of x and y in (6.13) and (6.14) will be exactly 
reversed. 



example 6.2 The Brachistochrone 




A famous problem in the calculus of variations is this: Given two points 1 and 
2, with 1 higher above the ground, in what shape should we build a frictionless 
roller coaster track so that a car released from point 1 will reach point 2 in 
the shortest possible time? This problem is called the brachistochrone problem, 
from the Greek words brachistos meaning “shortest” and chronos meaning 

I “time.” The geometry of the problem is sketched in Figure 6.4, where I have 

| 

taken point 1 as the origin and I have chosen to measure y vertically down. 

The time to travel from 1 to 2 is 

I f 2 ds 

1 time(I —> 2) = / — (6.17) 

i J 1 v 

I 

j where the speed at any height y is determined by conservation of energy to be 
j v = yigy. (Problem 6.8.) Because this gives v as a function of y , it is convenient 
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Figure 6.4 The brachistochrone problem is to find the shape of 
track on which a roller coaster released from point 1 will reach 
point 2 in the minimum possible time. 


to take y as our independent variable. That is, we shall write the unknown path 
as x = x (y). This means that the distance ds between neighboring points on the 
path has to be written as 


ds = V dx 1 + dy 2 = y/x'(y) 2 + 1 dy (6.18) 

where a prime now denotes differentiation with respect to y; that is, x'(y ) = 
dx/dy. Thus according to (6.17) the time of interest is 


time(l 



1 r 2 ^x'(y) 2 + 1 d 
V2 g Jo Vy 


(6.19) 


Equation (6.19) gives the integral whose minimum we have to find. It is of 
the standard form (6.14), except that the roles of x and y have been interchanged, 
with the integrand 

fix, x', y) = - t --, (6.20) | 

Vy [ 

To find the path that makes the time as small as possible, we have only to apply | 
the Euler-Lagrange equation (again with x and y interchanged) to this function, j 


d£ = ±df_ 

dx dy dx' 


( 6 . 21 ) 


The function of (6.20) is independent of x, so the derivative df/dx is zero, and 
(6.21) tells us simply that df/dx' is a constant. Evaluating this derivative (and 
squaring it for convenience) we conclude that 


x 


a 


y(l + x' 2 ) 


1 

const = — 
2 a 


( 6 . 22 ) 


where I have named the constant 1/2 a for future convenience. This equation is 
easily solved for x' to give 
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! whence 



This integral can be evaluated by the unlikely looking substitution 

y = a(l — cos#) 

which gives (as you should check) 


x — a J(1 — cos #) dO 
= a {6 — sin#) + const. 


(6.23) 


(6.24) 


(6.25) 


The two equations (6.25) and (6.24) are parametric equations for the required 
path, giving x and y as functions of the parameter #. We have chosen the initial 
point 1 to have x = y = 0, so we see from (6.24) that the initial value of # is 
zero. This in turn implies that the constant of integration in (6.25) is zero. Thus 
the final parametric equation for the path is 


x = a (9 — sin#) and y — < 2(1 — cos#) (6.26) 

with the constant a chosen so the curve passes through the given point (x 2 , y 2 ). 

The curve (6.26) is plotted in Figure 6.5. In that figure I have continued 
the curve (with dashes) beyond the point 2 to show that the curve that solves 
the brachistochrone problem happens to be a cycloid — the curve traced out by 
a point on the rim of a wheel of radius < 2 , rolling along the underside of the 
x axis (Problem 6.14). Another remarkable feature of this curve is this: If we 
release the cart from rest at point 2 and let it roll to the bottom of the curve 
(point 3 in the figure), the time to roll from 2 to 3 is the same whatever the 
position of 2, anywhere between 1 and 3. This means that the oscillations of a 
cart rolling back and forth on a cycloid-shaped track are exactly isochronous 
(period perfectly independent of amplitude), in contrast with the oscillations of 
a simple pendulum, which are only approximately isochronous, to the extent 

I 



Figure 6.5 The path for a roller coaster that gives the shortest 
time between the given points 1 and 2 is part of the cycloid with a 
vertex at 1 and passing through 2. The cycloid is the curve traced 
by a point on the rim of a wheel of radius a that rolls along the 
underside of the x axis. Point 3 is the lowest point on the curve. 
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that the amplitude is small. (See Problem 6.25.) The isochronous property of 
the cycloid was actually used in the design of some clocks, one of which can be 
seen in the Victoria and Albert Museum in London. 


Maximum and Minimum vs. Stationary 

You have probably noticed that in neither example of this section did I check that the 
curves that we found actually gave a minimum value to the integral of interest — that 
the straight line between two points actually makes the path length minimum , not a 
maximum or just stationary. The Euler-Lagrange equation guarantees only to give a 
path for which the original integral is stationary. The problem of deciding whether we 
have a minimum or maximum (or a stationary curve that is neither) is generally very 
difficult. In a few cases, it is easy to see which is the case. For instance, it really is 
obvious that a straight line gives the minimum distance between two points in a plane. 
In the case of the brachistochrone, it is not at all obvious that the path we found does 
yield a minimum time, though it is in fact true. 

To illustrate the variety of possibilities, consider the problem of finding the shortest 
path, or geodesic, between two points 1 and 2 on the surface of a globe. As you 
probably know, the answer is the great circle joining the two points. 4 Using the calculus 
of variations you can prove relatively easily that a great circle does indeed make the 
distance stationary: Using spherical polar coordinates, every point on the globe can 
be identified by the two angles 0 and 0. If you characterize a path as 0 = 0 (0) and set 
up an integral that gives the distance between 1 and 2 along this path, you can show 
that the Euler-Lagrange equation for 0(0) requires that the path follow a great circle. 
(See Problem 6.16 for details.) But you have to think a little carefully before deciding 
that this necessarily gives a minimum distance, since there are two different great- 
circle paths connecting any two points 1 and 2 on the globe: For simplicity consider 
two towns on the equator, Quito (near the Pacific coast of Ecuador) and Macapa (at 
the mouth of the Amazon on the Atlantic coast of Brazil). The “right” shortest path 
between these two is, of course, the great-circle path following the equator for about 
2000 miles across South America. But a second possibility, which satisfies the Euler- 
Lagrange equation just as well, is to head west around the equator from Quito, across 
the Pacific, the African continent, and the Atlantic, arriving in Macapa some 22,000 
miles later. You might guess that this path would be a maximum, but it is in fact neither 
maximum nor minimum: It is easy to construct nearby paths that are shorter, but it is 
also easy to find others that are longer. In other words, this second great-circle path 
gives neither a maximum nor a minimum. This second path is, of course, analogous 
to the horizontal point of inflection in elementary calculus. In this problem, luckily, 
it is obvious that the first path gives the true minimum. However, it should be clear 
that, in general, deciding what sort of stationary path the Euler-Lagrange equation 
has given us can be tricky. 


A great circle is the circle in which the globe intersects a plane through the globe’s center. 
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Fortunately for us, these questions are irrelevant for our purposes. We shall find that 
for the applications in mechanics all that matters is that we have a path which makes 
a certain integral stationary. It simply doesn’t matter whether it gives a maximum, 
minimum, or neither. 


6.4 More than Two Variables 


So far we have considered only problems with just two variables, the independent 
variable (usually x) and the dependent (usually y ). For most applications in mechanics, 
we shall find that there are several dependent variables, though fortunately still only 
one independent variable, which is usually the time t. For a simple example where 
there are two dependent variables, we can go back to the problem of the shortest path 
between two points. When we found the shortest path between two points 1 and 2, we 
assumed that the required path could be written in the form y = y (x). Reasonable as 
this seems, it is easy to think of paths that cannot be written in this way, such as the 
path shown in Figure 6.6. If we want to be perfectly sure we have found the shortest 
path among all possible paths, we must find a method that includes these. The way to 
do this is to write the path in parametric form as 

x = x(u) and y = y(u), (6.27) 

where u is any convenient variable in terms of which the curve can be parameterized 
(for instance, the distance along the path). The parametric form (6.27) includes all 
of the curves considered before. [If y = y(jc), just use x for the parameter u.] It also 
includes curves like that of Figure 6.6 and, in fact, all curves of interest. 5 

The length of a small segment of the path (6.27) is 

ds = y/dx 2 + dy 2 = yjx\u) 2 + y'(u) 2 du (6.28) 

where, as usual, a prime denotes differentiation with respect to the function’s argu¬ 
ment; that is, x\u ) = dx/du and y\u) = dy/du. Thus the total path length is 

L — I yfx'(u) 2 + y'(u) 2 du , (6.29) 

J U\ 

and our job is to find the two functions x(u) and y(u) for which this integral is 
minimum. 

This problem is more complicated than any we have considered before, because 


5 In case you are interested in mathematical niceties, I should say that in what follows I shall 
assume that all functions concerned are continuous and have continuous second derivatives. This 
assumption could be weakened a little, for example, by allowing some discontinuities in the 
derivatives. 
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Figure 6.6 This path between the two points 1 and 2 cannot be 
written as y = y(x) nor as i = x(y). It can be written in the 
parametric form (6.27). 


there are now two unknown functions x ( u ) and y(u). The general problem of this type 
is this: Given an integral of the form 


S 



f[x(u ), y(u ), x\u), y'(u ), u]du 


(6.30) 


between two fixed points [x(u{), y(wi)] and [x(u 2 ), y(u 2 )], find the path [ x(u ), y(u)] 
for which the integral S is stationary. The solution to this problem is very similar to 
the one-variable case, and I shall just sketch it, leaving you to fill in the details. The 
upshot is that with two dependent variables, we get two Euler-Lagrange equations. 
To prove this, we proceed very much as before. Let the correct path be given by 


x=x(u) and y = y(u). 


(6.31) 


and then consider a neighboring “wrong” path of the form 

x = x(u) + a£(u) and y — y(u) + /3r}(u) (6.32) 


(where c is the Greek letter “xi”). The requirement that the integral S be stationary for 
the right path (6.31) is equivalent to the requirement that the integral S(a, ft), taken 
along the wrong path (6.32), satisfy 


— = 0 and — = 0 (6.33) 

da dp 

when a = p = 0. These two conditions are the natural generalization of the condition 
(6.10) for the one-variable case. By an argument which exactly parallels that leading 
from (6.10) to (6.13), you can show that these two conditions are equivalent to the 
two Euler-Lagrange equations (see Problem 6.26): 


9/ d df 


3 y du 3 y' 


dx du dx' 


and 


(6.34) 
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These two equations determine a path for which the integral (6.30) is stationary, and, 
conversely, if the integral is stationary for some path, that path must satisfy these two 
equations. 


example 6.3 The Shortest Path between Two Points Again 

We can now solve completely the problem of the shortest path between two 
points. (That is, solve it including all possible paths, such as that in Figure 6 . 6 .) 
From (6.29), we see that for this problem the integrand / is 


f{x, x', y, y', u ) = V *' 2 + y a . 


(6.35) 


Since this is independent of x and y, the two derivatives df/dx and df/dy on 
the left sides in (6.34) are zero. Therefore, the two Euler-Lagrange equations 
imply simply that the two derivatives df/dx' and df/dy' are constants, 


9/ 

dx r 


x' 


\J x' 2 + y' 2 


= Ci 


and 


y 


9/ 


dy' yjx' 2 + y’ 2 


— C 2 . (6.36) 


If we divide the second equation by the first and recognize that y'/x' is just the 
derivative dy/dx , we conclude that 


dy 

dx 


y 

X' 


Cy 

C, 


= m, 


(6.37) 


say. It follows that the required path is a straight line, y = mx + b. It is interest¬ 
ing that this proof using a parametric equation is not only better than our previous 
j proof (in that the new proof includes all possible paths), it is also marginally 

i 

easier. 


The generalization of the Euler-Lagrange equation to an arbitrary number of 
dependent variables is straightforward, and doesn’t need to be spelled out in detail. 
Here I would just like to sketch the way the Euler-Lagrange equations will appear in 
the Lagrangian formulation of mechanics. 

The independent variable in Lagrangian mechanics is the time t. The dependent 
variables are the coordinates that specify the position, or “configuration,” of a system, 
and are usually denoted by q h q 2 , • • • ,q n • The number n of coordinates depends on the 
nature of the system. For a single particle moving unconstrained in three dimensions, 
n is 3, and the three coordinates q x , q 2 , q 3 could be just the three Cartesian coordinates 
x , y, z, or they might be the spherical polar coordinates r, 0, (p.Fov N particles moving 
freely in three dimensions, n is 3N and the coordinates q\,'",q n could be the 3 N 
Cartesian coordinates x h y h z h - • -, x N , y N , z N . For a double pendulum (two simple 
pendulums, with the second suspended from the bob of the first, as in Figure 6.7), 
there would be two coordinates q\,q 2 , which could be chosen to be the two angles 
shown in Figure 6.7. Because the coordinates q h • • •, q n can take on so many guises, 
they are often referred to as generalized coordinates. It is often useful to think of 
the n generalized coordinates as defining a point in an n -dimensional configuration 
space, each of whose points labels a unique position, or configuration, of the system. 
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Figure 6.7 A good choice of generalized coordinates to identify 
the position of a double pendulum is the pair of angles 0 X and 0 2 
between the pendulums and the vertical. 


The ultimate goal in most problems in Lagrangian mechanics is to find how the 
coordinates vary with time; that is, to find the n functions q x (t ), • • •, q n (t). One can 
regard these n functions as defining a path in the n-dimensional configuration space. 
This path is, of course, determined by Newton’s second law, but we shall find that it 
can, equivalently, be characterized as the path for which a certain integral is stationary. 
This means that it must satisfy the corresponding Euler-Lagrange equations (called 
just Lagrange equations in this context), and it turns out that these Lagrange equations 
are usually much easier to write down and use than Newton’s second law. In particular, 
unlike Newton’s second law, Lagrange’s equations take exactly the same simple form 
in all coordinate systems. 

The integral S whose stationary value determines the evolution of the mechanical 
system is called the action integral. Its integrand is called the Lagrangian L and 
depends on the n coordinates q x , q 2 , • • •, q n , their n time derivatives q h q 2 , • • •, q n 
and the time t , 


L=L{q x ,q x ,---,q n ,q n ,t). (6.38) 

Notice that since the independent variable is t, the derivatives of the coordinates q t 
are time derivatives and are denoted, as usual, with dots as q t . The requirement that 
the action integral 


rh 

S= &(qi, q h • • •, q n , q n , t) dt (6.39) 


be stationary implies n Euler-Lagrange equations 


3L d dL dL d dL , dL d dL 

— =--, -=-—, •••, and -=- 7 -. (6.40) 

dq x dt dq x dq 2 dt dq 2 3 q n dt 3 q n 

These n equations correspond precisely to the two Euler-Lagrange equations in (6.34) 
and are proved in exactly the same way. If these n equations are satisfied, then the 
action integral (6.39) is stationary; and if the action integral is stationary, then these n 
equations are satisfied. In the next chapter, you will see where these equations come 
from and how to use them. 
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Principal Definitions and Equations of Chapter 6 

The Euler-Lagrange Equation 


An integral of the form 


S = 


f f[y(x),y'(x),x]dx 
J X i 


[Eq. (6.4)] 


taken along a path y = y(x) is stationary with respect to variations of that path if and 
only if y(x) satisfies the Euler-Lagrange equation 


dy dx dy r 


= 0 . 


[Eq. (6.13)] 


Several Variables 

If there are n dependent variables in the original integral, there are n Euler-Lagrange 
equations. For instance, an integral of the form 


S = 


/*« 2 
J U\ 


f[x(u) 9 y(u ), x'(u), y'(u ), u]du, 


with two dependent variables [x (u) and y (m)], is stationary with respect to variations 
of x(u) and y(u) if and only if these two functions satisfy the two equations 


9/ 

dx 


d df 


du dx 1 


and V 


d 3 / 


dy du dy' 


[Eq. (6.34)] 


Problems for Chapter 6 _ 

Stars indicate the approximate level of difficulty, from easiest (*) to most difficult (+++). 


section 6.1 Two Examples 

6.1 ★ The shortest path between two points on a curved surface , such as the surface of a sphere, is 
called a geodesic. To find a geodesic, one has first to set up an integral that gives the length of a path on 
the surface in question. This will always be similar to the integral (6.2) but may be more complicated 
(depending on the nature of the surface) and may involve different coordinates than x and y. To illustrate 
this, use spherical polar coordinates (r, 0, 0) to show that the length of a path joining two points on a 
sphere of radius R is 


L = R j^ 2 rf+^erfe) 2 de (6.41) 


if (0 1? ff} and (0 2 , $ 2 ) specify the two points and we assume that the path is expressed as (j) = 4>(0). 
(You will find how to minimize this length in Problem 6.16.) 
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Figure 6.8 Problem 6.3 


6.2 ★ Do the same as in Problem 6.1 but find the length L of a path on a cylinder of radius R , using 
cylindrical polar coordinates (p, 0, z). Assume that the path is specified in the form 0 = 0(z). 

6.3 ★★ Consider a ray of light traveling in a vacuum from point P x to P 2 by way of the point Q on a 
plane mirror, as in Figure 6.8. Show that Fermat’s principle implies that, on the actual path followed, Q 
lies in the same vertical plane as P x and P 2 and obeys the law of reflection, that 0 l = 0 2 . [Hints: Let the 
mirror lie in the xz plane, and let P x lie on the y axis at (0, y b 0) and P 2 in the xy plane at (x 2 , y2>°)- 
Finally let Q = (x, 0, z). Calculate the time for the light to traverse the path P\QP 2 and show that it is 
minimum when Q has z = 0 and satisfies the law of reflection.] 

6.4 ★★ A ray of light travels from point P x in a medium of refractive index n ] to P 2 in a medium of 
index n 2 , by way of the point Q on the plane interface between the two media, as in Figure 6.9. Show 
that Fermat’s principle implies that, on the actual path followed, Q lies in the same vertical plane as P x 
and P 2 and obeys Snell’s law, that n x sin 0 X = n 2 sin0 2 - [Hints: Let the interface be the xz plane, and 
let P x lie on the y axis at (0, h h 0) and P 2 in the x, y plane at (x 2 , — h 2 , 0). Finally let Q = (x, 0, z). 
Calculate the time for the light to traverse the path P\QP 2 and show that it is minimum when Q has 
z = 0 and satisfies Snell’s law.] 

6.5 ★★ Fermat’s principle is often stated as “the travel time of a ray of light, moving from point A 
to B, is minimum along the actual path.” Strictly speaking it should say that the time is stationary , 
not minimum. In fact one can construct situations for which the time is maximum along the actual 
path. Here is one: Consider the concave, hemispherical mirror shown in Figure 6.10, with A and B 
at opposite ends of a diameter. Consider a ray of light traveling in a vacuum from A to B with one 
reflection at P, in the same vertical plane as A and B. According to the law of reflection, the actual 



Figure 6.9 Problem 6.4 
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^0 

Figure 6.10 Problem 6.5 


path goes via point P 0 at the bottom of the hemisphere (0 = 0). Find the time of travel along the path 
APB as a function of 0 and show that it is maximum at P — P 0 . 

6.6 ★★ In many problems in the calculus of variations, you need to know the length ds of a short segment 
of a curve on a surface, as in the expression (6.1). Make a table giving the appropriate expressions for 
ds in the following eight situations: (a) A curve given by y = y(x) in a plane, (b) same but x = x(y), 
(c) same but r = r (</>), (d) same but 0 = 0 (r); (e) curve given by 0 = 0 (z) on a cylinder of radius R , 
(f) same but z = z( 0 ); (g) curve given by 0 = 0 ( 0 ) on a sphere of radius R , (h) same but 0 = 0 ( 0 ). 

section 6.3 Applications of the Euler-Lagrange Equation 

6.7 ★ Consider a right circular cylinder of radius R centered on the z axis. Find the equation giving 0 
as a function of z for the geodesic (shortest path) on the cylinder between two points with cylindrical 
polar coordinates (/?, 0 1? z { ) and ( R , 0 2 , z 2 ). Describe the geodesic. Is it unique? By imagining the 
surface of the cylinder unwrapped and laid out flat, explain why the geodesic has the form it does. 

6.8 ★ Verify that the speed of the roller coaster car in Example 6.2 (page 222) is y/2gy. (Assume the 
wheels have negligible mass and neglect friction.) 

6.9 ★ Find the equation of the path joining the origin O to the point P( 1, 1) in the xy plane that makes 
the integral {y' 2 + yy 7 + y 2 ) dx stationary. 

6.10 ★ In general the integrand /(y, y 7 , jc) whose integral we wish to minimize depends on y, y 7 , and x. 
There is a considerable simplification if / happens to be independent of y, that is, / = /(y 7 , x). (This 
happened in both Examples 6.1 and 6.2, though in the latter the roles of x and y were interchanged.) 
Prove that when this happens, the Euler-Lagrange equation (6.13) reduces to the statement that 

df/dy' = const. (6.42) 

Since this is a first-order differential equation for y (x), while the Euler-Lagrange equation is generally 
second order, this is an important simplification and the result (6.42) is sometimes called a first integral 
of the Euler-Lagrange equation. In Lagrangian mechanics well see that this simplification arises when 
a component of momentum is conserved. 

* 2 /- 

<s/x-sj\ + y f2 dx is stationary. 

i 

6.12 ★★ Show that the path y = y(x) for which the integral f x^/l — y' 2 dx is stationary is a sinh 
function. 


6.11 ★★ Find and describe the path y = y(x) for which the integral 
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6.13 ** In relativity theory, velocities can be represented by points in a certain “rapidity space” in 
which the distance between two neighboring points is ds = [2/(1 — r 2 )]^/dr 2 + r 2 d0 2 , where r and 
0 are polar coordinates, and we consider just a two-dimensional space. (An expression like this for 
the distance in a non-Euclidean space is often called the metric of the space.) Use the Euler-Lagrange 
equation to show that the shortest distance from the origin to any other point is a straight line. 

6.14 ★★ (a) Prove that the brachistochrone curve (6.26) is indeed a cycloid, that is, the curve traced 
by a point on the circumference of a wheel of radius a rolling along the underside of the x axis, 
(b) Although the cycloid repeats itself indefinitely in a succession of loops, only one loop is relevant 
to the brachistochrone problem. Sketch a single loop for three different values of a (all with the same 
starting point 1) and convince yourself that for any point 2 (with positive coordinates x 2 , y 2 ) there is 
exactly one value of a for which the loop goes through the point 2. (c) To find the value of a for a given 
point x 2 , y 2 usually requires solution of a transcendental equation. Here are two cases where you can 
do it more simply: For x 2 = nb, y 2 = 2b and again for x 2 = 2nb, y 2 = 0 find the value of a for which 
the cycloid goes through the point 2 and find the corresponding minimum times. 

6.15 ★★ Consider again the brachistochrone problem of Example 6.2 (page 222) but suppose that the 
car is launched from point 1 with a fixed speed v 0 . Show that the path of minimum time to the fixed 
point 2 is still a cycloid, but with its cusp (the top point of the curve) a height v 2 /2g above point 1. 

6.16 ★★ Use the result (6.41) of Problem 6.1 to prove that the geodesic (shortest path) between two 
given points on a sphere is a great circle. [Hint: The integrand /(</>, 0', 0) in (6.41) is independent of 0, 
so the Euler-Lagrange equation reduces to 3//30' = c, a constant. This gives you 0' as a function of 0. 
You can avoid doing the final integral by the following trick: There is no loss of generality in choosing 
your z axis to pass through the point 1. Show that with this choice the constant c is necessarily zero, 
and describe the corresponding geodesics.] 

6.17 ★★ Find the geodesics on the cone whose equation in cylindrical polar coordinates is z = Xp. [Let 
the required curve have the form 0 = 0(p).] Check your result for the case that X —> 0. 

6.18 ★★ Show that the shortest path between two given points in a plane is a straight line, using plane 
polar coordinates. 

6.19 ★★ A surface of revolution is generated as follows: Two fixed points (xj, yj) and (x 2 , y 2 ) in the 
x, y plane are joined by a curve y = y(x). [Actually you’ll make life easier if you start out writing 
this as x = x(y).] The whole curve is now rotated about the x axis to generate a surface. Show that 
the curve for which the area of the surface is minimum has the form y = y 0 cosh[(x — x 0 ) /y 0 ], where 
x 0 and y 0 are constants. (This is often called the soap-bubble problem, since the resulting surface is 
usually the shape of a soap bubble held by two coaxial rings of radii y x and y 2 .) 

6.20 ★★ If you haven’t done it, take a look at Problem 6.10. Here is a second situation in which you 
can find a “first integral” of the Euler-Lagrange equation: Argue that if it happens that the integrand 
/(y, y r , x) does not depend explicitly on x, that is, / = /(y, y'), then 


df 

dx 


3 / , 3 / „ 

= —y H-y . 

3y 3 y' 


Use the Euler-Lagrange equation to replace 3//3y on the right, and hence show that 
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This gives you the first integral 

/ — /— = const. (6.43) 

dy r 

This can simplify several calculations. (See Problems 6.21 and 6.22 for examples.) In Lagrangian me¬ 
chanics, where the independent variable is the time t , the corresponding result is that if the Lagrangian 
function is independent of t, then energy is conserved. (See Section 7.8.) 

6.21 ★★ In Example 6.2 (page 222) we found the brachistochrone by exchanging the variables x and 
y. Here is a method that avoids that exchange: Write the time as in Equation (6.19) but using x as the 
variable of integration. Your integrand should have the form /(y, y', x) = >/(y /2 + 1 )/y. Since this is 
independent of x, you can invoke the “first integral” (6.43) of Problem 6.20. Show that this differential 
equation leads you to the same integral for x as in Equation (6.23) and hence to the same curve as 
before. 

6.22 ★★★ You are given a string of fixed length / with one end fastened at the origin O , and you are to 
place the string in the xy plane with its other end on the x axis in such a way as to enclose the maximum 
area between the string and the x axis. Show that the required shape is a semicircle. The area enclosed 
is of course f y dx, but show that you can rewrite this in the form f l {) fds , where 5 denotes the distance 
measured along the string from O , where / = yj 1 — y /2 , and y r denotes dyIds. Since / does not 
involve the independent variable s explicitly, you can exploit the “first integral” (6.43) of Problem 
6.20. 

6.23 ★★★ An aircraft whose airspeed is v Q has to fly from town O (at the origin) to town P, which 
is a distance D due east. There is a steady gentle wind shear, such that v wind = Vy x, where x and y 
are measured east and north respectively. Find the path, y = y(x), which the plane should follow to 
minimize its flight time, as follows: (a) Find the plane’s ground speed in terms of u 0 , V, 0 (the angle 
by which the plane heads to the north of east), and the plane’s position, (b) Write down the time of 
flight as an integral of the form fdx. Show that if we assume that y' and 0 both remain small (as 
is certainly reasonable if the wind speed is not too large), then the integrand / takes the approximate 
form / = (1 + \y' 2 )/(l + ky) (times an uninteresting constant) where k = V/v 0 . (c) Write down the 
Euler-Lagrange equation that determines the best path. To solve it, make the intelligent guess that 
y(x) = Xx(D — x), which clearly passes through the two towns. Show that it satisfies the Euler- 
Lagrange equation, provided X = (V4 + 2 k 2 D 2 — 2 )/(kD 2 ). How far north does this path take the 
plane, if D = 2000 miles, v 0 = 500 mph, and the wind shear is V =0.5 mph/mi? How much time does 
the plane save by following this path? [You’ll probably want to use a computer to do this integral.] 

6.24 ★★★ Consider a medium in which the refractive index n is inversely proportional to r 2 ; that is, 
n = a/r 2 , where r is the distance from the origin. Use Fermat’s principle, that the integral (6.3) is 
stationary, to find the path of a ray of light travelling in a plane containing the origin. [Hint: Use two- 
dimensional polar coordinates and write the path as 0 = 0(r). The Fermat integral should have the 
form //(0, 0', r) dr , where /(</>, 0', r) is actually independent of 0 . The Euler-Lagrange equation 
therefore reduces to 3// 30' = const. You can solve this for 0' and then integrate to give 0 as a function 
of r. Rewrite this to give r as a function of 0 and show that the resulting path is a circle through the 
origin. Discuss the progress of the light around the circle.] 

6.25 ★★★ Consider a single loop of the cycloid (6.26) with a fixed value of a, as shown in Figure 6.11. 
A car is released from rest at a point P 0 anywhere on the track between O and the lowest point P (that 
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Figure 6.11 Problem 6.25 


is, P Q has parameter 0 < 0 o < n). Show that the time for the cart to roll from P 0 to P is given by the 
integral 


time 


P) = 




1 — cos 0 
cos 0 o — cos 0 


dO 


and prove that this time is equal to n^/a/g. Since this is independent of the position of P Q , the cart 
takes the same time to roll from P Q to P, whether P 0 is at O, or anywhere between O and P, even 
infinitesimally close to P. Explain qualitatively how this surprising result can possibly be true. [Hint: 
To do the mathematics, you have to make some cunning changes of variables. One route is this: Write 
0 = 7 T — 2a and then use the relevant trig identities to replace the cosines of 0 by sines of a. Now 
substitute sin a — u and do the remaining integral.] 


section 6.4 More than Two Variables 

6.26 ★★ Give in detail the argument that leads from the stationary property of the integral (6.30) to the 
two Euler-Lagrange equations (6.34). 

6.27 ★★ Prove that the shortest path between two points in three dimensions is a straight line. Write 
the path in the parametric form 

x=x(u), y = y(u), and z = z(u) 

and then use the three Euler-Lagrange equations corresponding to (6.34). 



